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The notion of singular reduction modules, i.e., of singular modules of nonclassical (conditional) 
symmetry, of differential equations is introduced. It is shown that the derivation of nonclassical 
symmetries for differential equations can be favourably enhanced by an in-depth prior study of 
the associated singular modules of vector fields. The form of differential functions and differential 
equations possessing parameterized families of singular modules is described up to point trans- 
formations. Singular cases of finding reduction modules are related to lowering the order of the 
corresponding reduced equations. As examples, singular reduction modules of evolution equations 
and second-order quasi-linear equations are studied. Reductions of differential equations to alge- 
braic equations and first-order ordinary differential equations are considered in detail within the 
I ^ , framework proposed. 

1 Introduction 

^P" 1 ! The "nonclassical" approach to finding solutions of differential equations in closed form was 

' ■^ | proposed in [4J using the particular example of the (l+l)-dimensional heat equation in order to 

extend the range of applicability of symmetry methods. Since the end of the 1980s this method 
was applied to many particular differential equations modeling real- world phenomena, see, e.g., 
examples in [3] El \7\ [26] and reviews in [HEI]. Related objects, which are similar to subalgebras 
of Lie symmetry algebras, are named in the literature in different ways: nonclassical [TB], Q- 
conditional |9j, conditional |12j . partial [29] symmetries for short, or involutive families/modules 
\ of nonclassical/conditional symmetry operators |21[ IB2] in a more complete form. The main 

feature which is inherited by nonclassical symmetries from Lie symmetries is that they allow to 
construct ansatzes for the unknown function which reduce the corresponding differential equation 
to differential equations with a smaller number of independent variables \19\ \T7\ l29| [32] . This 
' feature relates nonclassical symmetries to the direct method by Clarkson and Kruskal [5j and 

the general ansatz method [Qj. In fact, however, properties of nonclassical symmetries are more 
closely connected with theories of differential constraints and formal compatibility of systems of 
differential equations |14[ \19\ [27] . This is why we mostly use the term "reduction modules" (of 
vector fields) instead of "involutive families of conditional symmetry operators" and say that an 
involutive module of vector fields reduces a differential equation if the equation is reduced by 
the corresponding ansatz. 

Involving the associated invariant surface condition in the conditional invariance criterion 
gives rise to a few significant complications of nonclassical symmetries in comparison with Lie 
symmetries. Given a differential equation C, elements of its different reduction modules form 
no objects of a nice algebraic or differential structure. Hence it is not possible to compose single 
reduction operators in reduction modules as this is done for the maximal Lie invariance algebra 
of C and its subalgebras, which consist of vector fields generating one-parametric Lie symmetry 
pseudogroups of C. Whereas the system of determining equations for Lie symmetries is linear, 
similar systems for reduction modules are nonlinear and should additionally be supplemented, in 
the course of considering modules of dimension greater than one, by the condition of involutivity, 
i.e., the closure of modules with respect to commutation of vector fields. Moreover, there is 
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no single system of determining equations even for reduction modules of a fixed dimension. 
Instead, the entire set of such modules is partitioned into subsets associated with systems of 
determining equations which are rather different from each other. Solving some of these systems 
may be equivalent to solving the initial equation, which gives no-go cases of looking for reduction 
modules. Such no-go cases were known for a number of particular (l+l)-dimensional evolution 
equations including the linear heat equation [8j [TUl E3 [30] , the Burgers equation [H |T7] , linear 
second-order evolution equations |23t [25] as well as for the entire classes of (l+l)-dimensional 
evolution equations [31] . multi-dimensional evolution equations [23] and even systems of such 
equations [28] . Note that in the course of the study of Lie symmetries a similar no-go situation 
arises for first-order ordinary differential equations [18} Section 2.5]. In fact, all the above no-go 
cases of reduction operators are occurrences of a no-go case common for evolution equations and 
one more no-go case specific for linear second-order evolution equations. The causes giving rise 
to the partition of the module set and to no-go cases for nonclassical symmetries had not been 
investigated in the literature until recently. It was not understood in what way results on no-go 
cases can be extended to reduction modules of other, non-evolution, equations. 

In [15] the partition of the set of reduction modules of a differential equation was related 
with lowering the order of this equation on the manifolds determined by the associated invariant 
surface conditions in the appropriate jet space. As a result, studying singular modules of vector 
fields which lower the order of the equation was included as the initial step in the procedure of 
finding nonclassical symmetries. In order to explicate the main ideas of the framework proposed, 
we considered only the case of single partial differential equations in two independent and one 
dependent variables and single reduction operators. The notion of singular reduction operators 
was introduced. The weak singularity co-order of a reduction operator Q were shown to be equal 
to the essential order of the corresponding reduced equation and the number of essential param- 
eters in the family of Q-invariant solutions. No-go assertions on singular reduction operators of 
(l+l)-dimensional evolution and wave equations were derived and then generalized to param- 
eterized families of vector fields which reduce partial differential equations in two independent 
variables to first-order ordinary differential equations. 

In the present paper we extend results of [15] to the case of a greater number of independent 
variables. After giving the basic notions and statements on nonclassical symmetries in Section [21 
in Section [3] we introduce the concepts of singular and meta-singular modules of vector fields 
for differential functions. Any meta-singular module of dimension greater than two proves to be 
necessarily involutive, in contrast to two-dimensional meta-singular modules. The main result 
of this section is Theorem [TJ which describes, up to point transformations, differential functions 
possessing meta-singular modules. The analogous notions of weakly singular and meta-singular 
modules for differential equations are introduced in Section 01 Theorem [21 which character- 
izes differential equations admitting weakly meta-singular modules, implies that instead of such 
modules it suffices to study meta-singular modules of the corresponding differential functions. A 
connection between the weak singularity co-order of reduction operators, the essential order of 
the corresponding reduced equations and, in the case of reduction to ordinary differential equa- 
tions, the number of parameters in the corresponding families of invariant solutions is established 
in Section [5l Revisiting results of [13] within the framework of singular reduction modules, in 
Section [6] we consider the specific case of reduction modules of dimension equal to the number 
of independent variables, which results in the reduction to algebraic equations. In Section [7] 
we reformulate and extend no-go results from [24] on modules reducing evolution equations to 
ordinary differential equations with time as the single independent variable. This motivates the 
consideration of reduction modules of singularity co-order one in Section [H Supposing that a 
differential equation £ admits an n-dimensional meta-singular module M of singularity co-order 
one, where n is the number of independent variables in £, we prove no-go assertions establishing 
connections between (n— l)-dimensional reduction modules of C contained in M and solutions 
of C. In particular, it is shown that the determining equations for such modules are reduced to 
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the initial equation C by a composition of a differential substitution and a hodograph transfor- 
mation. The final Section [9] is devoted to singular modules for quasi-linear second-order PDEs. 
Thus, elliptic equations possess no singular modules. Any evolution equation whose matrix 
of coefficients of second-order derivatives is non-degenerate possesses only singular modules as 
considered in Section [7] for general evolution equations. Generalized wave equations are much 
more complicated from this point of view. In particular, they may admit families of singular 
modules which have no interpretation in terms of meta-singular modules, which makes a further 
development of the framework of singular modules desirable. 

2 Reduction modules of differential equations 

In this section, based on [TTJ [121 ED [261 E2] , we briefly recall the required notions and results on 
nonclassical (conditional) symmetries of differential equations. After substantiating our choice 
with different arguments, we use the name "reduction modules" instead of "involutive families 
of nonclassical (conditional) symmetry operators" . 

Given a foliated space of n independent variables x = (x\, . . . ,x n ) and a single dependent 
variable u, consider a finite-dimensional involutive module Q of vector fields in this space, and 
suppose that the module dimension p of Q (over the ring of smooth functions of (x,u)) is 
not greater than n, < p ^ n. We additionally assume that the module Q satisfies the rank 
condition, i.e., for each fixed value of (x, u) the projection of Q to the space of x is p-dimensional. 
The attribute 'involutive' means that the commutator of any two vector fields from Q belongs 
to Q. It is obvious that any one-dimensional module is involutive. Therefore, in the case p = 1 
we can omit the attribute 'involutive' and talk only about modules. 

In what follows the indices i and j run from 1 to n, the index s runs from 1 to p, and 
we use the summation convention for repeated indices. Angular brackets (...) are used for 
denoting linear spans over the ring of smooth functions of (x, u). Subscripts of functions denote 
differentiation with respect to the corresponding variables, di = d/dxi and d u = d/du. Any 
function is considered as its zero-order derivative. All considerations are in the local setting. 

Suppose that the vector fields Q s = ^ Sl (x,u)di + n s (x,u)d u form a basis of Q, i.e., Q = 
(Ql, . . . , Q p )- Then the rank condition is equivalent to the equality rank(£ SJ ) = p. The condition 
that the commutator of any pair of basis elements belongs to Q, [Q S ,Q S '] £ Q, suffices for the 
module Q to be involutive. If the vector fields Q\, . . . , Q p form another basis of Q then there 
exists a nondegenerate p x p matrix-function (X ss (x,u)) such that Q s = X ss Q s /. 

The first-order differential function Q s [u] := n s (x,u) — £ st (x,u)ui is called the characteristic 
of the vector field Q s . In view of the Frobenius theorem, involutivity of Q is equivalent to the 
fact that the characteristic system Q[u] = of PDEs Q s [u] = (also called the invariant surface 
condition) has n + 1 — p functionally independent integrals u°(x, u), . . . , u n ~ p (x, u). Therefore, 
the general solution of this system can implicitly be represented in the form F(uj°, . . . , uj n ~ p ) = 0, 
where F is an arbitrary smooth function of its arguments. 

A differential function = @[z] of the dependent variables z = {z , . . . ,z m ) which in turn 
are functions of a tuple of independent variables y = (yi, . . . ,y n ) will be viewed as a smooth 
function of y and derivatives of z with respect to y. The order r = ordG of the differential 
function equals the maximal order of derivatives involved in 0. The differential function 
is defined as a function on a subset of the jet space J r = J r (y\z) with independent variables y 
and dependent variables z [IB] . 

Using another basis of Q gives just another representation of the characteristic system 
Q[u] = with the same set of solutions. This is why the characteristic system Q[u] = is 
associated rather with the module Q than with a fixed basis of Q. And vice versa, any family 
of n + 1 — p functionally independent functions of x and u is a complete set of integrals of the 
characteristic system of an involutive p-dimensional module. Therefore, there exists a one-to-one 
correspondence between the set of involutive p-dimensional modules and the set of families of 
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n + 1 — p functionally independent functions of x and u, which is factorized with respect to the 
corresponding equivalence. (We consider two families of the same number of functionally inde- 
pendent functions of the same arguments as equivalent if any function from one of the families 
is functionally dependent on functions from the other family.) 

A function u = f(x) is called invariant with respect to the involutive module Q (or, briefly 
Q-invariant) if it is a solution of the characteristic system Q[u] = 0. This notion is justified 
by the following facts. In view of the rank condition, we can choose a basis of Q which spans, 
over the underlying field, a p-dimensional (Abelian) Lie algebra g of vector fields in the space 
(x,u). The graph of each solution of the characteristic system Q[u] = is obviously invariant 
with respect to the p-parametric local transformation group generated by the algebra q. 

In view of the rank condition we can assume without loss of generality that ^ and 
F u o and resolve the equation F = with respect to co°: lo° = ^(co 1 , ■ ■ ■ ,Lu n ~ p ). This 
representation of the function u is called an ansatz corresponding to the module Q. 

Next, consider an rth order differential equation C of the form L(x, u^ r) ) = for the unknown 
function u of the independent variables x = (x%, . . . ,x n ). Here, u (r) denotes the set of all the 
derivatives of the function u with respect to x of order not greater than r, including u as the 
derivative of order zero. In the local approach the equation C can be viewed as an algebraic 
equation in the rth order jet space J r = J r (x\u) and is identified with the manifold of its 
solutions in J r : 

C = {(x, U( r )) G J r | L(x, ur r \) = 0}. 

We use the same symbol C for this manifold and write Q (r) for the manifold defined by the set 
of all the differential consequences of the characteristic system Q[u] =0 in J r , i.e., 

Q (r) = {(x,u {r) ) G J r \D^ ■■■D° n Q 8 [u] = 0, Oj G NU {0}, \a\: = a x + ■ ■ ■ + a n < r}, 

where D$ = d Xi +u ot j r $ i d Ua is the operator of total differentiation with respect to the variable Xi, 
a = (ax, ■ ■ ■ , o. n ) is an arbitrary multi-index, and 5i is the multi- index whose ith entry equals 1 
and whose other entries are zero. The variable u a of the jet space J r corresponds to the derivative 
d^u/dxl 1 . . .dx" n , and Ui = ug v 

Definition 1. The differential equation C is called conditionally invariant with respect to 
the involutive module Q if the relation V^ r - ) L(x,u^ r) )\cnQ ir) = holds for any V G Q. This 
relation is called the conditional invariance criterion. Then Q is called an involutive module of 
conditional symmetry (or Q-conditional symmetry, or nonclassical symmetry, etc.) operators of 
the equation C. 

In this definition, Vr r \ denotes the standard rth prolongation of a vector field V = Ci x i u )di + 
r](x,u)d u [M [22]: V (r) = V + Z 0< \a\<rV a du a , where rf = ■ ■ ■ D%»V[u] + Cu a+&1 and 
V[u] = rj — t^Ui. It suffices to check the conditional invariance criterion for V running through 
a basis of Q. What basis is chosen for representing the characteristic system Q[u] =0 and 
checking the conditional invariance criterion is not essential |12t [32] . 

The equation C is conditionally invariant with respect to the module Q if and only if an ansatz 
constructed with this module reduces C to a differential equation with n—p independent variables 
(resp. to an algebraic equation if n = p) [32J. Thus, we will briefly call involutive modules of 
conditional symmetry operators of C reduction modules of C. The set of p-dimensional reduction 
modules of the equation C will be denoted by 1Z P (£). 

An alternative approach to conditional invariance is to demand that the joint system of C 
and Q (r ) is formally compatible in the sense of the absence of nontrivial differential conse- 
quences [19} 121] . If the conditional invariance criterion is not satisfied but nevertheless the 
equation L has Q-invariant solutions then one talks about weak invariance of the equation C 
with respect to the module Q [201 [2TT [27] . 
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There are reduction modules related to classical Lie symmetries. Let be a p-dimensional 
Lie invariance algebra of the equation C, whose basis operators Qi, . . . , Q p satisfy the condition 
rank(£ Si ) = rank(^ sl , r/ s ) ( = p' ^ p). Then the span of Q\, . . . , Q p over the ring of smooth 
functions of (x, u) is a p'-dimensional involutive module which belongs to TZ P {C). Modules of 
this kind are called Lie reduction modules. Other reduction modules are called non-Lie. 

A useful alternative formulation of the conditional invariance criterion is as follows (cf. [32J). 

Lemma 1. Given an rth order differential equation C: L[u] = 0, a p- dimensional (0 < p < n) 
involutive module Q satisfying the rank condition and differential functions L[u] and X[u] ^ 
of an order not greater than r such that L\q m = \L\q,., the module Q is a reduction module 

of £ if and only if it is a reduction module of the equation C: L[u] = 0. 

The classification of reduction modules can be notably enhanced and simplified by involving 
Lie symmetry and equivalence transformations of (classes of) differential equations. By 5UT P 
we denote the set of p-dimensional modules of vector fields in the space of (x,u). Any point 
transformation of (x, u) induces a one-to-one mapping of 9Jl p into itself via push- forward of 
vector fields. Namely, the transformation g: Xj = X l (x,u), u = U(x,u) generates the mapping 
g* : S0t p TV such that for any Q G TV and V G Q the vector field V = u)di + rj(x, u)d u is 
mapped to the vector field g^V = + rjd%, where u) = VX l (x,u), fj(x,u) = VU(x,u). 

Given a group G of point transformations in the space of (x, u), the modules Q and Q (of the 
same dimension) are called equivalent with respect to G if there exists some g G G such that 
Q = g*Q- 

Lemma 2. Suppose that Q G 1Z P {£), a point transformation g maps an equation C to an 
equation C and the image g*Q satisfies the rank condition. Then g*Q G 1Z P (£). 

Corollary 1. Let G be the point symmetry group of an equation C Then the equivalence of p- 
dimensional modules of vector fields with respect to the group G generates an equivalence relation 
in TZ p {C). 

Next, we consider a class C\s of equations Cg: L(x,um,9) = parameterized by the 
parameter-functions 9 = 6(x,ur r \). Here L is a fixed function of x, ur r \ and 9. By 9 we denote 
the tuple of arbitrary (parametric) differential functions 9(x,u^) = (9 1 (x,u^.- ) ), . . . ,9 k (x,u^)) 
traversing the set S of solutions of a system. This system consists of differential equations 
S(x, ur r \, 0( g )(x, U( r ))) = and differential inequalities E(x, ur r \,6r q \(x, ur r \)) ^ (> 0, < 0, . . . ) 
on 9, where x and Ur r \ play the role of independent variables and 9r q \ stands for the set of all 
the derivatives of 9 of order not greater than q. Henceforth we call the functions 9 arbitrary 
elements. We write G~ for the point transformation group preserving the form of the equations 
from C\s- 

By P we denote the set of all pairs consisting of an equation Cg from C\s and a module Q 
from TZ p (£g). It follows from Lemma [2] above that the action of transformations from the 
equivalence group G~ on £\$ and {TV \ 9 G S} induces an equivalence relation on P |26j. 

Definition 2. Let 9,9' G S, Q G K p (£ e ), Q' G H p {C e >). The pairs (C g ,Q) and (C e >,Q') 
are called G~ -equivalent if there exists g G G~ such that g transforms the equation Cq to the 
equation Cq>, and Q' = g*Q. 

We will interpret the classification of reduction operators with respect to G~ as the clas- 
sification in P with respect to this equivalence relation, a problem which can be investigated 
similarly to the usual group classification in classes of differential equations. Namely, at first 
we construct the reduction modules that are defined for all values of 9. Then we classify, with 
respect to G~, the values of 9 for which the equation Cg admits additional reduction modules. 
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3 Singular modules of vector fields for differential functions 



Let L = L[u] be a differential function of order ordL = r (i.e., a smooth function of the 
independent variables x = (x\, . . . , x n ) and derivatives of u with respect to x up to order r) and 
let Q be a p-dimensional (0 < p < n) involutive module which is generated by the vector fields 
Qs = £, st {x.,u)di + 7] s (x,u)d u defined in the space of (x,u) and satisfying the rank condition 
rank(£ sl ) = p. 

Definition 3. The module Q is called singular for the differential function L if there exists a 
differential function L = L[u] of an order less than r such that L\q,. = L\q. Otherwise Q 
is called a regular module for the differential function L. If the minimal order of differential 
functions whose restrictions on Q (r) coincide with L\q,. equals k (k < r) then the module Q 
is said to be of singularity co- order k for the differential function L. The module Q is called 
ultra- singular for the differential function L if L\q.. = 0. 

It will be advantageous to define the singularity co-order of ultra-singular modules and the 
order of identically vanishing differential functions to be —1. If a module is regular for the 
differential function L we define its singularity co-order to be r = ordL. The singularity co- 
order of the module Q for the differential function L will be denoted by scol Q- 

The case p = n is special. Given an n-dimensional involutive module Q which is generated 
by vector fields satisfying the rank condition, for any rth order differential function L = L[u] 
there exists a zero-order differential function L = L[u] such that L\q,. = L\q, . . This is why in 
this case it is natural to assume that the module Q is singular for L only if it is ultra-singular 
for L, i.e., L\ Q(r) = 0. 

Next we show how to algorithmically construct a function L satisfying the conditions of 
Definition [3l First note that without loss of generality we may suppose that rank(£ ss ) = p and 
consider the basis Q s = d s + ^ SL d L + fi s d u , where the index i runs from p + 1 to n and the matrices 
and (fj s ) are the products of the matrix (£ ss ) _1 by the matrices (£ st ) and (77 s ), respectively. 
Then any derivative of u of order not greater than r can be expressed, on the manifold Q(r), via 
derivatives of u with respect to x p+ i, . . . , x n only and the coefficients £ st and fj s . For example, 
for the first- and second-order derivatives we have 

u SL = fjf - ifu t > + fj s u u t - C'u L ,u L - i sl 'u u ,, 

Uss , = f) s s , - i s j Ui + - Cu L )(f' - £ sV u t ! ' 

Substituting the expressions for the derivatives into L, we obtain a differential function L de- 
pending only on x, u and derivatives of u with respect to x p+ \, . . . , x n . We will call L a 
differential function associated with L on the manifold Q( r j. The module Q is singular for the 
differential function L if and only if the order of L is less than r. The co-order of singularity 
of Q equals the order of L. The module Q is ultra-singular if and only if L = 0. This shows that 
testing for an involutive module satisfying the rank condition to be singular for a differential 
function is indeed feasible in an entirely algorithmic way and can easily be included in existing 
programs for symbolic calculations of symmetries. 

Proposition 1. Suppose that L = L[u] is a differential function of u = u{x), x = (xi, . . . , x n ), 
and let Q be an involutive module of vector fields defined in the space of (x,u), which is of 
dimension less than n and satisfies the rank condition. Then scol Q ^ sco^ Q for any involutive 
submodule Q of Q. In particular, the module Q is singular for L if it contains a submodule 
singular for L. 
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Proof. If Q is an involutive submodule of Q then it necessarily satisfies the rank condition and 
Qir) Q Q(r)i where r = ordL. If the differential function L coincides with a differential function L 
on the manifold Q(r), the same is true on the manifold Q (r) . Therefore, sco^ Q ^ scol Q. □ 

Definition 4. A (p+l)-dimensional module M is called meta-singular for the differential func- 
tion L if any p-dimensional involutive submodule of M which satisfies the rank condition is 
singular for L and the module M contains a family SCR = {Q*} °f such submodules parame- 
terized by an arbitrary function = $(x,u) of all independent and dependent variables. The 
singularity co-order of the meta-singular module M is the maximum of the singularity co-orders 
of its involutive p-dimensional submodules satisfying the rank condition. 

Here and in what follows the parameterization by an arbitrary function is understood as a 
parameterization modulo functions with a smaller number of arguments. The singularity co- 
order of the family 59? = {Q*} is also defined as the maximum of the singularity co-orders of 
its elements. 

The condition that the module M is involutive is not explicitly included in Definition [H 
By this definition, the meta-singular module M should only contain a family of p-dimensional 
involutive submodules parameterized by an arbitrary function of all variables. At the same time, 
in the case p ^ 2 this is equivalent to the fact that the module M is involutive. 

Proposition 2. A (p+1)- dimensional module M, where p ^ 2, contains a family of p- 
dimensional involutive submodules parameterized by an arbitrary function of all independent 
and dependent variables if and only if the module M is involutive. 

Proof. Suppose that the module M contains a family 9JT of p-dimensional involutive submodules 
parameterized by an arbitrary function of all variables modulo functions with a smaller number 
of arguments. Then we can choose a basis {Qq, ■ ■ ■ , Qp} of M in such a way that the vector 
fields Qi, . . . , Q p generate an involutive submodule of M from 9Jt and, moreover, commute. By 
a change of variables (x,u) these vector fields are reduced to shift operators, Q s = d s . Up to 
combining with Q s and multiplying by a nonvanishing function, the vector field Qq can be chosen 
in the form Qq = £ 0t (x, u)d L + n°(x, u)d u , where one of the coefficients £ 0t , t = p + 1, ■ ■ ■ , n, or 77 
is equal to 1. The entire set of p-dimensional submodules of M is partitioned into the subsets 

S = {Q S = (Q s + s Qo, s = l,...,p}}, 

S s = {Q Ss = (Qo, Qs-i, Q S ' + s ' Q s , s' = s + 1, . . . ,p)}, 

where = (<? ,..., P ), s = (0 S+1 , ■ ■ ■ , P ), and all 0's run through the set of smooth functions 
of By construction, if Q® s is an involutive module then the chosen basis elements com- 

mute. This implies that the components of s satisfy the equations Qq0 s =0, . . . , Q S -\0 S = 0, 
s' = s + 1, . . . ,p, and can hence be expressed via functions of at most n + l — s arguments. There- 
fore, the module M contains a family 9Jt of p-dimensional involutive submodules parameterized 
by an arbitrary function of all variables if and only if such a family is a subset of Sq. 

If a module Q e from So is involutive, by the Frobenius theorem the overdetermined system 
of the first-order linear partial differential equations Q s & + s Qo& = with respect to the 
unknown function $ = <£(x, u) has a solution with Qo& 7^ 0, i.e., the parameter-functions s can 

1 This can be seen as follows: Suppose that R = {Ri, . . . , R p ) is a p-dimensional submodule of M and write 
Rs = T,l=o X3k Qk, where X ak are smooth functions of (x,u). Set A := (A sfe ) and A ' : = (^t^Zp- We 
have rank A = p. If rank A' = p, by elementary row operations we may generate a basis of R that identifies it 
as an element of So- Otherwise we may transform A into a form with A = 1 and all other entries of the first 
column and first row vanishing. Applying the same reasoning as before to the resulting lower right submatrix 
we obtain either a basis of type Si or we may again simplify the first column and first row of this submatrix as 
before. Thus the claim follows by induction. 
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be represented in the form 9 s = —& s /Qo&. As the basis elements Q s + 9 s Qo of the module Q 9 
should commute, we have 

[Q s + 9 s Qo, Qs' + # s Qo] = (#! + s Qq9 s — 9 S S , — 9 s Qo9 s )Qo + 9 s Qo, s — 9 s Qo, s ' 

<f>s>Q ,s® - $sQo,s>$ „ $s n $ s ' n _ n 

" (QoW 2 Q ° + Q^ Qo ' s ' ~ Q^ Q °' S ~ ' 

where Q , s = [Q s , Q ] = ^ + rf s d u . 

Suppose that not all the vector fields Qq :S are zero. Let vector fields Qi, Qp form a 
basis of the module (Qo,s)- Then the vector fields Qo, Q±, . . . , Qp are linearly independent, the 
vector fields Qo,s are represented as ^f=i ^ ss Qs for some smooth functions X ss = X ss (x,u) and 
the above commutation relations imply the following system for the function <E>: 

$s'Qo,s$ ~ $sQo,s>$ = 0, \ S ' S <S>s = ^ SS <S>s>, 8 = 1,..., p. 

As some of the coefficients X ss are necessarily nonzero, in this case the function $ runs at 
most through the solution set of a system of first-order linear partial differential equations and, 
therefore, the number of its arguments is in fact less than n + 1. This contradicts the existence 
of a family of p-dimensional involutive submodules of M parameterized by an arbitrary function 
of all variables modulo functions with a smaller number of arguments. 

As a result, all the vector fields Qo,s = [Qs, Qo] are zero, i.e., the module M is involutive. 

Conversely, if the module M is involutive, we choose a basis which consists of commuting 
vector fields Qo, . . . , Q p . Each submodule Q® = (Q s — (Q s ^)/(Qo^)Qo) is involutive and of 
dimension p. Here $ = &(x, u) runs through the set of smooth functions of (x, u) with Qo& ^ 0, 
i.e., the complement of the solution set of the equation Qo& = 0. The general solution of 
this equation is parameterized by a single function of n arguments. The modules and 
associated with the functions $ and <£ coincide if and only if (Qq$)Q s Q = (Qs^)Qo^- We will 
call such functions <E> and $ equivalent. Thus, the set traversed by <3? should additionally be 
factorized with respect to this equivalence relation. A necessary and sufficient condition for the 
equivalence of functions and <I> is that there exists a smooth function F of n — p + 1 arguments 
for which <3> = Ffa 1 , . . . , uj n ~ p , <J>), where oj 1 , uj n ~ p form a complete set of functionally 
independent solutions of the system Qqoj = 0, Q s u = 0. As the number of arguments of F is 
less then n + 1, the factorization does not affect the degree of arbitrariness of <I>. □ 

Remark 1. Families of involutive submodules belonging to the set So can be parameterized in 
different ways. Thus, in the above proof these submodules are parameterized via the represen- 
tation of the coefficients 9 s in the form 9 s = —Q s &/Qo&, and the set traversed by <3? = <J>(x, u) 
should be factorized with respect to an equivalence relation. Alternatively, we can choose, e.g., 
any of the 9 s instead of as a parameterized function. After VP 1 = 9 1 is chosen as such a 
function, the condition \Q\ + 9 1 Qo,Q2 + 9 2 Qo] = implies the first-order linear partial dif- 
ferential equation Q\9 2 + 9 1 Qq9 2 = Q29 1 + 9 2 Qo9 1 with respect to 9 2 , whose general solution 
is parameterized by an arbitrary function ^ 2 of n arguments. In the same way, the condition 
[Qi + 9 % Qo, Qs + 9 3 Qo] = 0, i = 1, 2, gives the system of two first-order linear partial differential 
equations Qi9 3 + 9 l Qo9 3 = Q^9 l + 9 3 Qq9 % . This system is associated with the module generated 
by the vector fields Qi + 9 l Qo + (Qffi + 9 3 Qo9 i )dg3, which is involutive in view of the equation 
for 9 2 . It follows from the Frobenius theorem that the general solution of the system is parame- 
terized by an arbitrary function ^ 3 of n — 1 arguments. Iterating the procedure for each 9 s , we 
precisely parameterize involutive modules from the submodule set 5*0 by the tuple {^f 1 , ■ ■ . , ^ p ), 
where runs through the set of smooth functions of n + 2 — s arguments. 

Corollary 2. Any (p+l)- dimensional meta- singular module M for a differential function L, 
where p ^ 2, is involutive. 
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Proposition 3. A (p+1)- dimensional module M , where p ^ 2, is meta- singular for a differential 
function L if and only if it contains, up to point transformations in the space of(x, u), a family 9Jt 
of p- dimensional involutive submodules singular for L of the form 

Q* = {d s -($ s /$ u )d u ,s = l,...,p) 

parameterized by an arbitrary smooth function = &(x,u) of all independent and dependent 
variables with $ u ^0. The singularity co- order of the family 9Jt for the differential function L 
coincides with that of the entire module M , sco^ 9JT = sco^ M . 

Proof. Suppose that the module M is meta-singular for a differential function L. As the mod- 
ule M is involutive in view of Corollary [2j we can choose a basis of M consisting of commuting 
vector fields Qo, Q p such that the vector fields Q\, Q p satisfy the rank condition 
rank(£ s *) = p. A change of the variables (x,u) reduces these vector fields to shift operators, 
Q s = d s and Qo = d u . A submodule of the form Q e = (Q s + 9 s Qo, s = 1, . . . ,p) of M is 
involutive if and only if the coefficients 9 s can be represented in the form 9 s = — <& s /& u for 
some function $ = <&(x,u) with <Ev 7^ 0, cf. the proof of Proposition [2j Thus, the family 
9Jt = {Q® = (d s — (Q s /Q u )d u , s = 1, . . . ,p)} parameterized by the arbitrary function $ with 
& u 7^ is of the required form. 

Conversely, let the module M contain, up to point transformations in the space of (x,u), 
such a family Wl = {Q*} and let sco^OJt = k < r = ordL. Hence we have that sco^Q* ^ k 
for any allowed value of the parameter-function <I>. In the initial coordinates, the basis elements 
of a submodule Q* take the form Q s — (Q s &)/(Qo^)Qo, where Qo, Q p are commuting 
vector fields. It suffices to prove that any p-dimensional involutive submodule P of M, which 
satisfies the rank condition and does not belong to the family {Q*}, is singular for L. Up to 
a permutation of the vector fields Qi, . . . , Q p , a basis of P consists of the vector fields Qo 
and Q s i + s Qi, s' = 2, . . . ,p, where the coefficients 6 s are smooth functions of (x,u). For 
convenience, by a point transformation of the variables (x, u) we reduce the vector fields Qo, Q2, 
. . . , Q p and Q\ to the shifts operators d\, 82, ■ ■ ■ , d p and d u , respectively. As the submodule P 
is involutive, the coefficients 9 s possess the representation 9 s = —ty s i/ty u , where ^ is a smooth 
function of x-i, ■ ■ ■ , x n and u with ^> u ^ 0. Consider the family of involutive modules of the form 

Q £ = (Qo + eQi, Qs' + 9 S ' £ Q U s' = 2,...,p) 

parameterized by a constant e running through a neighborhood of zero. Here the coefficients 
9 s 6 are obtained via replacing the argument u in 9 s by u — sx\, 9 s £ = 9 s (x2, . . . ,x n ,u — ex\). 
For each nonzero value of e the module Q e belongs to the family {Q*}. This is obvious after 
the transition from the chosen basis to the basis {(Qi + e^Qo, Q s ' — £~ 1 9 s £ Qo}- Therefore, 
any module Q e with e 7^ is singular for the differential function L. Consider the differential 
function L £ which is associated with L on the manifold Q^ r) via the exclusion of the derivatives u a 
with a\ + • • • + cx p > from L using the equations u\ = e, u s i = 9 s 6 and their differential 
consequences. The function L £ is smooth in the totality of the parameter e, the variables x and 
derivatives of u with respect to . . . , x n . As ordL £ ^ k for any e / 0, the same statement 

is true for e = through continuity. This means that the submodule P = Q £ | e= o is singular 
for L, and sco^ P ^ k = sco^ 9JT. 

It is obvious that sco^M ^ sco^iDT. As any p-dimensional involutive submodule Q of M 
satisfies the inequality sco lQ ^ sco l 9R, we obtain that sco ;,¥ = sco l □ 

The case of two-dimensional meta-singular modules is special. As any one-dimensional mod- 
ule of vector fields is involutive, the fact that a two-dimensional module is meta-singular for a 
differential function does not imply that this module is involutive. More specifically, the con- 
clusion of Proposition [5] is not true if p = 1. This is why the reduced form of submodules of a 
two-dimensional meta-singular module depends on whether this module is involutive or not. 
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Proposition 4. A two-dimensional module M is meta-singular for a differential function L if 
and only if it contains, up to point transformations in the space of (x,u), a family 9JT of one- 
dimensional submodules singular for L with basis vector fields in a reduced form parameterized 
by an arbitrary smooth function 9 = 6(x,u). The reduced form is 



if the module M is involutive or not involutive, respectively. Here £ l = £ % (x, u), i = 3, . . . , n are 
fixed smooth functions. The singularity co-order of the family 9Jt for the differential function L 
coincides with that of the entire module M , sco^ 9K = sco^ M. 

In order to make the further consideration for p = 1 consistent with the case of p 2, the 
parameter-function 9 can be represented in the form 9 = — $i/$ u , where $ = $>(x,u) is an 
arbitrary smooth function with <I> U 7^ 0. 

Proof. If the module M is involutive, the proof is similar to that of Proposition [2j We therefore 
only consider the case when the module M is not involutive. 

Let the module M be meta-singular for a differential function L. We choose a basis {Qo, Qi} 
of M such that the vector field Qi satisfies the rank condition rank(^ 1 *) = 1 and reduce the 
vector field Qo by a change of the variables (x, u) to the shift operator with respect to u, 
Qo = d u . Up to permutation of the variables x\, . . . , x p , we can assume that £ n ^ 0. Then 
we replace the basis element Qi by (£ n ) _1 (Qi — f^Qo) i n order to set r/ 1 = and £ n = 1. As 
the module M is not involutive, the commutator [Qo, Q%] = £* 2 £?2 + • • • + £,}pd n does not vanish. 
Hence we can assume up to permutation of the variables X2, ■ ■ ■ , x p that 7^ 0. The change 
of variables x s = x s and u = £ 12 (x,u) with the simultaneous replacement of Qo by (£„ 2 ) _1 Qo 
reduces the basis elements of M to the form Qo = d u and Q\ = d\ + u&2 + £ 3 (?3 + • • • + £, n d n . 
Then the family SPT = {{Qi + 9Qo)} of one-dimensional submodules of M singular for L, where 
the parameter 9 = 9(x, u) runs through the set of smooth functions of all independent and 
dependent variables, is of the required form. 

Conversely, let the module M contain, up to point transformations in the space of (x,u), 
such a family 9JT with sco^OK = k < r = ordL. Hence we have that scol(Qi + 9Qo) ^ k for 
any value of the parameter- function 9. After returning to the initial coordinates, it suffices to 
prove that the submodule (Qo) of M is singular for L if the vector field Qo also satisfies the 
rank condition in these coordinates. For convenience we reduce Qo by a change of coordinates 
to the shift operator d\. Consider the family of modules of the form Q e = (Qo + e(Qi — £ n Qo)) 
parameterized by a constant e running though a neighborhood of zero. For each nonzero value 
of £ the module Q £ belongs to the family 9Jt as in this case we have Q £ = (Qi + (e -1 — £ n )Qo)- 
Therefore, any module Q e with e 7^ is singular for the differential function L. Consider the 
differential function L £ which is associated with L on the manifold Q^ r) via the exclusion of the 
derivatives u a with ct\ > from L using the equation u\ = e(r/ 11 — £ 12 U2 — • • • — £ u n ) and 
its differential consequences. The function L £ is smooth in the totality of the parameter e, the 
variables x and derivatives of u with respect to X2, ■ ■ ■ , x n . As the order of L £ is not greater 
than k for any nonzero e, the same statement is true for e = by continuity. This means that 
the submodule (Qo) = Q e | e =o is singular for L, and scol(Qo) ^ k = scol Wl. 

The set of one-dimensional submodules of M is exhausted by Qo and the elements of the 
family Wl. Hence sco^ 9Jt = sco^ M. □ 

Theorem 1. An rth order differential function L with one dependent and n independent vari- 
ables possesses a kth co-order meta-singular (p+l)- dimensional module of vector fields if and 
only if it can be represented, up to point transformations, in the form 



d 1 + 98 u or dt+ud2 + C 



a 



L — L{x, £l r ,k,p)i 



(2) 
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where Q r> k,p = (w Q , |a| ^ r, a p+ \ + ••• + «„ ^ k) , L Wa ^ for some u a with a p+ \ + • ■ ■ + a n = k, 

and uo a = u a or, only for the case p = 1, m a = D^ 2 ■ ■ ■ D% n (D\ + uD 2 + £ 3 D 3 -\ + t; n D n ) ai u 

for some fixed smooth functions £ l = u), i = 3, . . . , n. 

Proof. Suppose that a differential function L possesses a kth co-order (p+l)-dimensional meta- 
singular module M. Up to combining basis elements and change of variables, a basis of the 
module M consists of either the vector fields Q s = d s and Qq = d u or, if p = 1 and the 
module M is not involutive, the vector fields Q\ = d\ + ud 2 + £ 3 <93 + • • • + C n d n and Qq = d u , 
cf. Propositions [3] and 01 Although the form of the initial differential function L will also be 
transformed by the change of variables, for simplicity we will continue to use the old notations 
for all new values. 

We choose the family 9K = {Q* = {Q s — ($ s /& u )d u , s = 1, . . . of p-dimensional involutive 
submodules of M which are parameterized by an arbitrary function <J> = <E>(x, u). Then we fix an 
arbitrary point z° = (x°,u® r) ) € J r and consider the p-dimensional involutive modules from SOT 
for which z° £ Q* )5 where Q* } denotes the manifold determined by the module in the 
jet space J r . This condition implies that the values of the derivatives of <J? with respect to 
only x\, x n at the point (x°,u°), which contain differentiation with respect to some x s , 
are expressed via u^ r) and values of derivatives of <1? in (x°,u°), containing differentiation with 
respect to u. The latter values are not constrained. For instance, if the module M is involutive, 
we get 

<f> s (x°,u ) = -u°$ u (x°, U ), 

d> si (x°, u°) = -u Mx°,u ) - u°<Z> su (x°,u°) - u° s $ iu (x , u°) - vP s u^ uu {x°, u°), .... 

We introduce the new coordinates {xj, Lv a , \a\ ^ r} in J r instead of the standard ones 
{xi,u a , \a\ ^ r}. If the module M is involutive, this change of coordinates is just a re-notation 
of variables in order to guarantee consistency with the special case of non-involutive modules 
for p = 1. In the latter case, this is a valid change of coordinates since the Jacobian matrix 
{du a / du a i) is nondegenerate: it is a triangular matrix with all diagonal entries equal to 1 if we 
implement the following order of multi-indices: 

a -< (3 & |a| < |ft V (\a\ = |ft A«i< ft) V (\a\ = [ft A a\ = ft A a 2 < ft) V • • • . 

Denote by L the differential function obtained from L by the procedure of excluding, on the 
manifold Qf r) , the derivatives of u involving differentiations with respect to x s (see (flj). For a 
multi-index a = (ati, . . . , a n ) we set a = [a\, . . . , a p ) and d = (a p +i, . . . , a n ), i.e., a = (d, a). 
The symbol denotes the tuple of p zeros. As is a feth co-order singular module for L, the 
function L does not depend on the derivatives UfQ ^, \a\ = k + 1, . . . , r. We use this condition 
step- by-step, starting from the greatest value of |d| and re- writing the derivatives in the new 
coordinates of J r and in terms of L. In the course of this procedure, we take into account the 
equality uip = ipP[u] satisfied on the manifold Qf r ^ for each multi-index (5 with [ft ^ r. Here the 

differential function ^ = ^[u] is defined by the equality ^ = ■ ■ ■ Dt(Qf) Pl • • • {Qpf p u, 

and hence it is of |/3|th order and possesses the representation 

^ = mQt) 131 ■ ■ ■ (Qtf p u)u m + $>[u] 

with some differential function ipP = ip^[u] of order less than [ft. Therefore, for each a with 
a = and \a\ ^ r the chain rule implies that 

L Ua (z°)= Y, L^(z°)^ Ua (z ). (3) 

13: |^|<r,|/3|>|d| 



11 



Thus, in the new coordinates for each a with a = and |d| = r the equation L Ua (z°) = 
can be written in the form L Wa (z°) = 0. Indeed, in this case we have that ipu a = 1 if /3 = a and 
4>u a = otherwise. This completes the first step. 

In the second step we fix a value of a with a = and |d| = r — 1. As L Wfl (z°) = if $ = and 
|/3 1 = r, the summation multi- index in ([3]) with the fixed a can be assumed to run through the 
set B\ = {j3 | |/3| ^ 1, |/3| = r — 1}. The derivative is equal to 1, d u Qfu and for j3 = (6, a), 
j3 = (S s ,a) and all other values of /3 from Bi, respectively. Here 5 S is the p-tuple with the sth 
entry equal to 1 and the other entries equal to 0. Therefore, the equation L u ^ (z°) = implies 
that 

L "{6,&)( Z °) + L ^(S s ,&)( z0 )( d ^Qf U )\(x,u)=(x»,u^ = °- 

Note that d u Qfu = —(& s /<fr u ) u . We split with respect to the value & su (x°,u°) as it is uncon- 
strained. Thereby we arrive at the equations L w fl (z°) = and L u ^ a <i)( z0 ) = 0- 

Iterating this procedure, before the /xth step, \x € {3, ... ,r — k}, we derive the equations 
Lup(z ) = 0, where the multi-index /3 runs through values for which r — ^ + 2^|/3|^r and 
|/3 1 ^ r— |/3 1 . Then for each fixed value of a with a = and \a\ = r — the summation multi- 
index in © can be assumed to run through the set = {/3 | |/3| ^ \x — 1, |/3| = r — /i + 1}. For 
/? € 5^ the derivative ^equals d u (Qf ■ ■ ■ {Qpf p u if $ = a and is zero otherwise. Therefore, 
the equation t> U() (z ) = implies the condition 

£ ) ft • • • (Q?)NU= ( *v°) = o. 

|^KM-l,/3=A 

The values (d^ 1 ■ ■ ■ dp p & u )(x° , u°), |/3| ^ \i — 1, are unconstrained. Hence by splitting with 
respect to them or, equivalently, by splitting with respect to (d u (Qf Y 1 ■ ■ ■ (Q* f p u) L , , , u0 y 
|/3| ^ \i — 1, we obtain the equations L^^z ) = 0, \j3\ = r — [A + l and |/3| ^ /x — 1. 

Finally, after the (r — /c)th step we derive the system L UJa (z°) = 0, |d| > k and \a\ ^ r, which 
implies the condition ([2]). 

Conversely, let an rth order differential function L be of the form ([2]) (after a point transfor- 
mation). For an arbitrary smooth function $ = <&(x,u) with <& u ^ 0, we consider the involutive 
module generated by either the vector fields Qf = d s — (Q s /Q u )d u in the general case or 
the vector field Qf = d\ + u&2 + ^ 3 c?3 + • • • + £ n d n — (& s /§ u )d u in the special case with p = 1. 
Using L and Q*, we construct the differential function L* = L(x,Q r ^, p ), where 

n rAp = [u a = D%? ■ --D^(Qfr ■ ■ ■ (QtT P u, \a\ ^ k, \a\ ^ r). 

By construction, ordL* ^ k for all values of the parameter-function $ with <£ u ^ 0. Moreover, 
ord L* = k for almost all values of this parameter- function except those which satisfy a system 
of differential equations. As 

L\Qf = L\q<S> , 

where $ runs through the set of smooth functions of (x, u) with nonvanishing derivative with 
respect to u, the family 93T = {Q*} is a A;th co-order singular family of p-dimensional involutive 
modules for the differential function L in the new coordinates. We return to the old coordi- 
nates. In view of Proposition [3] if p ^ 2 or Proposition [4] if p = 1, the module of vector fields 
which contains the family 9K is a fcth co-order meta-singular (p+l)-dimensional module for the 
differential function L. □ 

Excluding the special case of two-dimensional non-involutive meta-singular modules, the 
result presented in Theorem [1] can be formulated in the following way: A differential function 
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with one dependent and n independent variables admits a kih. co-order meta-singular (p+1)- 
dimensional involutive module M if and only if it can be reduced by a point transformation of 
the variables to a differential function in which the total differentiation with respect to n — p 
fixed independent variables in each derivative of the dependent variable is of order not greater 
than k. 

Corollary 3. Any differential function with one dependent and n independent variables (not 
identically vanishing) admits no meta-singular (p+1)- dimensional (0 < p < n) module of sin- 
gularity co-order —1 (i.e., whose p- dimensional involutive submodules are ultra- singular for this 
differential function). 

Remark 2. Clearly a feth co-order meta-singular (p+l)-dimensional module M may contain 
p-dimensional involutive modules whose singularity co-orders are less than the singularity co- 
order of the entire module M. Consider a family Wl = {Q* = (Q s — (& S /& U )Q° , s = 1, . . . ,p)} 
of p-dimensional involutive submodules of M which are parameterized by an arbitrary function 
<3? = <J>(x, u) and assume that 50? is singular for a differential function L, and sco^ 9JT = k. Here 
the vector fields Q s and Q° are assumed to be reduced to the form presented in the beginning of 
the proof of Theorem [TJ Then the values of <3? for which sco^ < k are solutions of the system 

L Wa (x, n rAp ) {d u (Qfr ■ ■ ■ (Qt) ap u) =o, \&\ = k, 

\a\^r—k 

where L and £l r ,k,p are defined in Theorem [T] and its proof, respectively. In other words, the 
regular values of <f> associated with the submodules of the maximal singularity co-order k in SOT 
satisfy, for some a with \a\ = k, the inequality 

£ L Ua (x,n rAp )(d u (Qf r • • • (QtT p u) + o. 

\a\^r—k 

4 Singular modules of vector fields for differential equations 

An involutive module Q satisfying the rank condition is called ( strongly ) singular for a differ- 
ential equation C if it is singular for the differential function L[u] constituting the left hand side 
of the canonical representation L[u] =0 of the equation C. Usually we will omit the attribute 
"strongly" . 

As left hand sides of differential equations are defined up to multipliers which are nonvan- 
ishing differential functions, the conditions from Definition [3] can be weakened when considering 
differential equations. 

Definition 5. A p-dimensional (0 < p < n) involutive module Q which satisfies the rank 
condition is called weakly singular for the differential equation C: L[u] = if there exists a 
differential function L = L[u] of an order less than r and a nonvanishing differential function 
A = X[u] of an order not greater than r such that L\q.. = (AL)|q, Otherwise we call Q 
a weakly regular module for the differential equation C. If the minimal order of differential 
functions whose restrictions on Q (r) coincide, up to nonvanishing functional multipliers, with 
L|g (r) is equal to k (k < r) then the module Q is said to be weakly singular of co-order k for 
the differential equation C. 

An involutive module Q is considered to be weakly ultra-singular for the differential equation 
C: L[u] = if it is strongly ultra-singular for C. As in the case of strong regularity, weakly regular 
modules for the differential equation C are defined to have weak singularity co-order r = ordL. 
We write wsco£ Q for the weak singularity co-order of the module Q for the equation C. 

Strong singularity implies weak singularity and consequently weak regularity implies strong 
regularity. The weak singularity co-order is always less or equal and may be strictly less than 
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the strong singularity co-order. Thus strongly regular vector fields may be singular in the weak 
sense. As a simple example, consider the equation 3:2^111+^1^222 = e" 33 (7/3+14). It possesses the 
two-dimensional singular module {81,82) whose strong and weak singularity co-orders equal 2 
and 1, respectively. The same module (81,82) is strongly regular and is of weak singularity 
co-order 1 for the equation x%u\\ + X1U22 = e n33 (i/3 + u). 

Let L be the differential function associated with L on the manifold Q (r) by excluding, in view 
of equations defining Q( V ), those derivatives of it which contain differentiations with respect to x s . 
We suppose that L is of maximal rank in a derivative u a of the highest order k appearing in this 
differential function, i.e., L u ^ 7^ for some a with |d| = A: on the solution manifold of the equa- 
tion L = (see the notation in the proof of Theorem [1]). It then follows that the weak singularity 
co-order of Q for the equation C: L = equals the order k of L and, consequently, the strong 
singularity co-order of Q. It follows that in this case there is an entirely algorithmic procedure 
for testing whether an involutive module is weakly singular for a partial differential equation. 

The following assertion is proven in the same way as Proposition [TJ 

Proposition 5. Suppose that £ is a differential equation with respect to the unknown func- 
tion u of n independent variables x and let Q be an involutive module of vector fields defined 
in the space of (x,u), which is of dimension less than n and satisfies the rank condition. Then 
wsco£ Q ^ wsco£ Q for any involutive submodule Q of Q. In particular, the module Q is weakly 
singular for C if it contains a submodule weakly singular for C. 

The notion of meta-singular modules for differential equations is defined similarly to the case 
of differential functions. Thus, we have the following definition of modules which are meta- 
singular in the weak sense. 

Definition 6. A (p+l)-dimensional module M is called meta-singular for the differential equa- 
tion L in the weak sense if any p-dimensional involutive submodule of M which satisfies the 
rank condition is weakly singular for L and the module M contains a family 9JT = {Q®} of 
such submodules parameterized by an arbitrary function <3? = $(a;,ii) of all independent and 
dependent variables. The weak singularity co-order of the meta-singular module M coincides 
with the maximum of the weak singularity co-orders of its involutive submodules satisfying the 
rank condition. 

Theorem 2. An rth order differential equation C: L[u] = of maximal rank with one dependent 
and n independent variables possesses a kth co-order weakly meta-singular (p+l)- dimensional 
module of vector fields if and only if the differential function L can be represented, up to point 
transformations of variables, in the form 

L = A[u]L(x,n rAp ), (4) 

where A is a nonvanishing differential function of order not greater than r, L is a smooth function 
of x and £l r fi yV = (ui a , \ot\ ^ r, a p+ i + - • -+ot n ^ kj , L ula ^ for some uj a with a p+ \ + - ■ -+a n = k, 

and uj a = u a or, only for the case p = 1, u a = D^ 2 ■ ■ ■ D" n (Di + uD 2 + ( 3 D 3 H h £ n D n ) ai u 

for some fixed smooth functions £ l = ^ l {x,u), i = 3, . . . ,n. The value of k should be minimal 
among all possible representations of the differential function L of the form (dD . 

Proof. We will freely use the notations and definitions from the proof of Theorem [TJ 

To begin with, assume that the differential equation C: L[u] = is of maximal rank and 
admits a kth co-order weakly meta-singular p+l-dimensional module of vector fields. Up to 
point transformations and changes of module basis, it suffices to consider only the family 9Jt = 
{Q* = (Q s — {&s/&u)Q°, s = l,... ,p)} of p-dimensional involutive submodules of M which are 
parameterized by an arbitrary function $ = <J?(x, u). Here Q° is reduced to the shift operator 8 U 
and the vector fields Q s take either the form Q s = 8 S or, if p = 1 and the module M is not 
involutive, the form Qi = d\ + u&2 + ^83 + • • • + £, n 8 n . 
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Given any point z° = (:r , -u^ ) € C C J r we choose the modules from 9Jt = {Q } for which 
•2° € Q*j- It then follows that the values of the derivatives of & with respect to only x\, . . . , x n 
at the point (x°,u°), which contain differentiation with respect to some x s , can be expressed 
via u^ r) and values of derivatives of <1? in (x°,u°), containing differentiation with respect to u. 
These latter values are then unconstrained. 

We next derive from L the differential function L by excluding, on the manifold Qf r) , deriva- 
tives of u involving differentiations with respect to x s (see ([I])). The fact that for any function $ 
with $ u the module is at most of fcth co-order weak singularity for C leads to the 
condition L u .~ (zo) = 0, |d| = fc + 1, . . . , r. As in the proof of Theorem [1] we now iteratively 
employ this condition, starting out with the greatest value r of |d| and re- writing the derivatives 
in the new coordinates {xi,co a , \a\ ^ r} of J r and in terms of L. Thereby we obtain 

Luj a (z°) = 0, \a\ > k, \a\ < r, 

which is satisfied for any z° £ C Applying the Hadamard lemma to each of these equations and 
then simultaneously integrating them, we obtain (cf. the proof of Theorem 1 in [32]). 

Conversely, suppose that the rth order differential function L is of the form (after a point 
transformation). For an arbitrary smooth function <I> = <£(x,u) with & u 7^ 0, we consider the 
involutive module Q* generated by either the vector fields Qf = d s — (& S /Q u )d u in the general 
case or the vector field Qf = d\ + u&z + £ 3 $3 + • • • + £, n d n — ($ S /$ U )5 U in the special case with 
p = 1. Using L and Q*, we construct the differential function = L(x, & r kp), where 

n r , ktP = (co a = D^X ■ --D^iQfr ■ ■ ■ (Qt) ap u, \&\ ^ k, \a\ ^ r). 

By construction, the order of the differential function L* is not greater than k for all allowed 
values of the parameter- function <E>. Moreover, ordL* = k for almost all values of this parameter- 
function except those which satisfy a system of differential equations. As 

LU* = ALU* , 

where $ runs through the set of smooth functions of (x, u) with nonvanishing derivative with 
respect to u, the family 9JI = {Q*} is a kth co-order weakly singular family of p-dimensional 
involutive modules for the differential equation C in the new coordinates. We return to the 
old coordinates. In view of Proposition [3] if p ^ 2 or Proposition H] if p = 1, the module M of 
vector fields which is spanned by the family 9Jt is a kth co-order meta-singular (p+l)-dimensional 
module for the differential function L/A. Therefore, this module is a fcth co-order weakly meta- 
singular (j>+l)-dimensional module for the differential equation C □ 

Corollary 4. A differential equation C: L[u] = of maximal rank with one dependent and 
n independent variables possesses a kth co-order weakly singular p- dimensional module of vector 
fields (0 < p < n) if and only if this module is kth co-order strongly singular for £ (possibly 
in a representation differing from L[u] = by multiplication by a nonvanishing differential 
function of u). 

5 Singularity of reduction modules and parametric 
families of solutions 

Definition 7. An involutive module Q is called a singular reduction module of a differential 
equation C if Q is both a reduction module of C and a weakly singular module of C 

The following assertions are obtained as direct extensions of the corresponding results for the 
case n = 2 [15] . 
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Proposition 6. Let Q be a p- dimensional non-ultra- singular reduction module (0 < p < n) of 
an equation £. Then the weak singularity co-order of Q for £ coincides with the essential order 
of a corresponding reduced differential equation. 

Proof. By a point transformation we may achieve the situation that in the new variables the 
module Q has a basis {Q s = d s , s = 1, . . . ,p}. (Again we use the same notation for old and 
new variables). Then an ansatz constructed with Q is u = (f(co), where ip = (p(uj) is the 
new unknown function, ui = (ui, . . . , u; n _ p ) and U\ = x p +i, . . . , w n _ p = x n are the invariant 
independent variables. The manifold Q (r) is defined by the system u a = 0, where the multi-index 
a = {ot\, . . . , a n ) satisfies the conditions a\ + ■ ■ ■ + a p > 0, \a\ := a\ + • • • + a n ^ r = ordL. 

As Q £ 1Z P (£), there are differential functions A = X[ip] and L = L[ip] of an order not greater 
than r such that L\ u=ip ^ = XL (cf. [32]). The function A does not vanish and may depend on 
the variables x s as parameters. Modulo equivalence generated by nonvanishing multipliers we 
may assume that the function L is of minimal order f. Then the reduced equation £: L = has 
essential order f. 

Now since wsco£ Q = k it follows that there exists a strictly /cth order differential function 
L = L[u] and a nonvanishing differential function A = X[u] of an order not greater than r, which 
depend at most on x and derivatives of u with respect to Xp+i, . . . , x n such that Z/|q^ — XL. 

Suppose first that f would be less than k. In this case we can use A new = X\ u ^ v and 
L new = L\ u ^ v in the definition of weak singularity to arrive at the contradiction wsco£ Q ^ 
ordZ new = f < k. Thus f ^ k. (Here "y z" means that the value y should be substituted 
instead of the value z.) 

On the other hand, if f > k we conclude that XL = (XL)\ u=(p ^, where the variables x s play 
the role of parameters. Fixing a value x° s of x s for each s, we obtain the representation 



L 



u=<p(uj), x 3 



However, as ord L\ u =<f>(w), x =x° K k < f, this representation contradicts the condition that f is 
the essential order of the reduced equation C. 

We conclude that f = k. The inverse change of variables preserves the claimed property. □ 

The properties of ultra-singular modules of vector fields as reduction modules are obvious. 

Proposition 7.1) Any p-dimensional module Q of vector fields (0 < p ^ n) which is ultra- 
singular for a differential equation C is a reduction module of this equation. An ansatz con- 
structed with Q reduces £ to the identity. Therefore, the family of Q -invariant solutions of £ is 
parameterized by an arbitrary function of n — p Q-invariant variables. 

2) If Q is a p-dimensional module of vector fields and the family of Q-invariant solutions 
of £ is parameterized by an arbitrary function of n — p Q-invariant variables then Q is an 
ultra- singular module of vector fields for £. 

Given a partial differential equation £, consider a reduction module Q for £ of dimension n— 1, 
i.e., p = n — 1. Such a module reduces the equation £ to an ordinary differential equation whose 
essential order, in view of Proposition [6j is equal to wsco£ Q. This allows us to relate wsco£ Q 
with the maximal number of parameters in families of Q-invariant solutions of £. An exception 
occurs for the ultra-singular modules which reduce £ to identities. In this case we define the 
order of the reduced equation to be —1. 

Proposition 8. Suppose that Q is an (n— 1)- dimensional module of vector fields and a differen- 
tial function of minimal order, associated with the differential function L[u] on the manifold <2 {r) 
up to a nonvanishing multiplier (r = ordL), is of maximal rank in the highest order derivative 
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of u appearing in this differential function. Then any two of the following properties imply the 
third one. 

1) Q is a reduction module of the equation C: L = 0. 

2) The weak singularity co-order of Q for C equals I, where ^ Z ^ r. 

3) The equation C possesses an l-parametric family of Q-invariant solutions, and any Q- 
invariant solution of C belongs to this family. 

Proof. If Q is an (n— l)-dimensional non-ultra-singular reduction module of an equation C then 
the weak singularity co-order of Q for C equals the maximal number Nc.q of essential continuous 
parameters in families of Q-invariant solutions of C. Indeed, the weak singularity co-order 
of Q for C coincides with the essential order f of the reduced ordinary differential equation C 
associated with Q, f = wsco£ Q. The maximal number of essential continuous parameters in 
solutions of L also equals f. The substitution of these solutions into the corresponding ansatz 
leads to parametric families of Q-invariant solutions of C, and all Q-invariant solutions of C are 
obtained in this way. Therefore, Ncq = f. 

In view of the condition of maximal rank, the equation C can be written in normal form and 
hence has an f-parametric general solution which contains all solutions of C. Substituting it 
into the corresponding ansatz, this solution gives an f-parametric family of Q-invariant solutions 
of C. There are no other Q-invariant solutions of C. Therefore, conditions 2 and 3 are equivalent 
if condition 1 holds. 

For any (n— l)-dimensional involutive module Q of vector fields we have Nc,q ^ k, where 
by k we denote wsco£ Q. Let us prove that Q is a reduction module of C if Nc q = k. 

Point transformations of the variables do not change the claimed property. We use the 
variables and notations from the proof of Proposition [6J Consider the differential function 
L[(p] = L\ u=ip r u y It depends on the variables x s as parameters, and ordL = k. Due to the 
condition of maximal rank, we can resolve the equation L = with respect to the highest order 
derivative <p( k ': ip^ = R[<p], where ordi? < k. 

If R Xa ^ for some s, splitting with respect to x s in the equation L = results in an 
ordinary differential equation R[<p] =0 of an order lower than k. Any Q-invariant solution of C 
has the form u = (f(uj), where the function (p satisfies, in particular, the equation R[(p] = 0. 
This contradicts the condition that Nc,q = k. 

Therefore, R Xa = 0, i.e., the equation ip( k > = R[(p] is a reduced equation obtained from the 
equation C by the substitution of the ansatz u = f{oS) constructed with the module Q. □ 

6 Reduction to algebraic equations 

As remarked in Section [31 the case when the dimension of modules of vector fields coincides 
with the number of independent variables, p = n, is special for the singularity of modules for 
differential functions. This case is also special for reduction of differential equations. In contrast 
to reduction modules of lower dimensions, n-dimensional reduction modules of any differential 
equation C with n independent variables reduce this equation to algebraic equations instead of 
differential equations with a smaller number of independent variables. Moreover, this is the only 
case when both regular and singular reduction modules can be studied in a uniform way. 

Let Q be an involutive module of dimension p = n which satisfies the rank condition. Then 
we can choose the basis of Q consisting of the vector fields Q s = d s +n s (x, u)d u . As the module Q 
is involutive, the basis elements Q s should commute, i.e., the coefficients n s = r] s (x,u) satisfy 
the system 

Vj + v'ti = V S s ' + V s Vu- (5) 

The Frobenius theorem implies that the system Q s & = & s + r ? s ^ > w = with respect to the 
function $ = <&(x,u) has a single functionally independent solution, which is not constant. In 
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other words, the coefficients rj s can be represented in the form 77 s = — <f> s /<3? u for some smooth 
function = $(2;, u) with & u 7^ 0. An implicit ansatz constructed for u with the module Q is 
<&(x,ii) = <p, where <p is the new unknown function. It is nullary since the module dimension p 
of Q equals the number n of independent variables x. 

Suppose that Q is a reduction module of an rth order differential equation C: L[u] = 0. 
All the derivatives of u with respect to x from order 1 up to order r are expressed, on the 
manifold Q( r ), via the variables x and u: 

u Q = h a (x, u) := (d 1 + r, l d u ) ai ■■■{8 n + v n d u ) an u, 1 < |a| < r. 

As the vector fields Q s commute, the above representation for derivatives of u is well defined 
since it does not depend on the ordering of operators in the right hand side expression. Using 
these expressions for excluding the appropriate derivatives u a from L, we get a differential 
function L = L[u] of order zero, i.e., a function of (x,u). Varying r/ s , it can also be interpreted 
as a differential function with the independent variables x and u and the dependent variables rj s . 
Therefore, in view of Lemma [T] the fact that Q is a reduction module of C is equivalent to the 
condition Q S L\^ =0 = 0, which holds in the zero-order jet space. This condition jointly with 
equations ([5]) gives the complete system of determining equations for the coefficients rj s . Under 
the representation r] s = — <3? S /<I> U it can be interpreted as an equation for the function where 
L* = L\ v s = _$ s /$ u is considered as a differential function with the independent variables x and u 
and the dependent variable <3>. 

In order to handle the condition Q S L\^ =Q = 0, we can solve the equation L = with respect to 
a variable (either one of the x's or u) and excluding this variable from the equation Q S L = using 
the obtained expression. As the functions L and r] s (resp. $) depend on the same arguments, 
the above procedure does not result in a usual differential equation. 

This is why we use another, more convenient, way, which involves the Hadamard lemma. 
Within the local approach, it suffices to split the further consideration into two cases: either the 
function L is of maximal rank or it is a constant. 

In the first case, for each s the condition Q S L\^ =0 = is equivalent, in view of the Hadamard 
lemma, to the equality Q S L = X S L for some smooth function A s = X s (x,u). After cross differ- 
entiation for each pair (s,s') we obtain 

Q S Q S ,L - Q S ,Q S L = Q S {\ S '~L) - Q S >(\ S L) = (Q S X S> - Q S 'X S )L = 0, 

i.e., Q S X S = Q S 'X S . Therefore, there exists a smooth function A = A(x,u) such that A s = Q S A- 
The system of the equations Q S L = X s L then implies in view of the representation ij s = —& s /<& u 
that = AC($) for some smooth function £ of one argument, where A = e is a nonvanishing 
function of (x,u). 

In the second case, the condition Q S L = is satisfied by each (x, u) as a differential con- 
sequence of the current supposition that the function L is a constant. Setting A = 1 and 
C{&) = = const, we derive the same equation L* = A£(<&) with the nonvanishing multi- 
plier A. 

In fact, this equation is precisely the condition of reduction of the equation C by the ansatz 
$(x,u) = cp associated with the module Q. Indeed, if the function u = u(x) is implicitly 
defined by this ansatz then the derivatives of u (of nonzero orders) are found from differential 
consequences of the equations u s = — $ s /$ u . Therefore, the substitution of the ansatz to C 
leads to the equation L \$r x ,u)=tp = which is equivalent, in view of the reduction condition, to 
the algebraic equation ((ip) = with respect to the constant (p. 

It is obvious that the module Q is ultra-singular for C if and only if the function £ identically 
vanishes. 

Summing up the above considerations results in the following assertion: 
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Proposition 9. A module Q is an n- dimensional reduction module of a differential equation C: 
L[u] = with n independent variables x if and only if this module is spanned by the vector fields 
d s — (&s/®u)du! s = 1, . . . ,n, where the function <1> = $>(x, u) satisfies the equation = A£(<&), 
A is a nonvanishing function of (x,u) and the differential function = L*[$] is obtained 
from L[u] by the exclusion of the derivatives of u using differential consequences of the equations 
u s = —§ S /Q U . The ansatz §(x,u) = ip reduces the equation C to the algebraic equation Cif) = 
with respect to the constant ip. 

Theorem 3. Up to the equivalence of solution families, for any differential equation C with 
respect to a single unknown function of n independent variables there exists a one-to-one corre- 
spondence between one-parametric families of its solutions and its n-dimensional ultra- singular 
reduction modules. Namely, each module of the above kind corresponds to the family of solutions 
which are invariant with respect to this module. The problems of the construction of all one- 
parametric solution families of the equation C and the exhaustive description of its n-dimensional 
ultra- singular reduction modules are completely equivalent. 

Proof. Let Q be an n-dimensional ultra-singular reduction module of the equation C. It follows 
from Proposition [9] that the ansatz <£(x,u) = tp constructed with the module Q reduces the 
equation C to the identity. In other words, for each value of the constant <p this ansatz implicitly 
defines a solution of C 

Conversely, suppose that T = {u = f(x,x)} is a family of solutions of C parameterized 
by the single constant parameter x. As this parameter is essential and, therefore, the deriva- 
tive f x is nonzero, we can express x from the equality u = f(x, x). As a result, we obtain that 
x = $(x, u) for some function $ = u) with <£ u ^ 0. Consider the module Q = (Qi, . . . , Q n ), 
where Q s = d s + r/ s d u and the coefficients r/ s = rj s (x,u) are defined by rj s = — $ S /<I> U . It is an 
n-dimensional involutive module and Q[u] = for any element of T. The ansatz u = f(x,(p), 
where ip is the new unknown (miliary) function, is associated with Q and reduces C to the 
identity. This means that Q is an ultra-singular reduction module for C 

One-parametric families T = {u = f(x, x)} and J 7 = {u = f(x, x)} are defined to be equiva- 
lent if they consist of the same functions and differ only by parameterizations, i.e., if there exists 
a function £ = C( x ) such that C, x ^ and f(x,C(x)) = f(x, x). This is true if and only if the 
functions $ = <&(x,u) and $ = $(x,u) associated with the families J- and J- ', respectively, are 
functionally dependent. More precisely, $ = C(^)- As ^ u ^u 7^ 0, the functional dependence of $ 
and <I> is equivalent to the equality $ s /$ u = Therefore, equivalent one-parametric fami- 

lies of solutions correspond to the same ultra-singular reduction module Q of C and, conversely, 
any two one-parametric families of Q-invariant solutions are equivalent. □ 

Corollary 5. The system of the determining equations for coefficients of ultra- singular reduc- 
tion modules of C, which consists of equations and the equation L = 0, is reduced by the 
composition of the nonlocal substitution r/ s = — § S /Q U , where $ is a function of (x,u), and the 
hodograph transformation 

the new independent variables: Xi = Xj, x = <3?, 
the new dependent variable: u = u 

to the initial equation C in the function u = u(t, x, x) with x playing the role of a parameter. 

Note that the reduction of differential equations to algebraic ones using nonclassical symme- 
tries was considered in [T3]. An assertion similar to Theorem was obtained therein. Using the 
notion of singular reduction operators makes Theorem [3] more precise. 
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7 Motivating example: evolution equations 



We investigate re-dimensional singular reduction modules of (l+n)-dimensional evolution equa- 
tions of the general form 



for the unknown function u depending on the variables t = xq and x = (x\, . . . , x n ). (For 
convenience, in this section we set the number of independent variables to n + 1 instead of n 
and additionally single out the variable xq.) Here it^) denotes the set of all derivatives of the 
functions u with respect to the space variables x of order not greater than r, including u as 
derivative of order zero. Also, ut = du/dt and r = max{|a| | H Ua ^ 0} ^ 2, i.e., we assume 
the order of the equations under consideration to be not less than two. We fix an arbitrary 
equation £ of the form ©. 

To reduce the equation C to an ordinary differential equation, we should use its re-dimensional 
reduction module. In general, such a module Q is spanned by re vector fields of the form 



with rank(r s , = n. (We use the same convention for the ranges of indices as in the whole 
paper, i.e., the indices i and j run from 1 to re and the index s runs from 1 to p noting that here 
p coincides with re). In contrast to the specific case re = 1 |15j . for general values of n we can 
completely describe singular n-dimensional modules of the equation C only after more precisely 
knowing the form of H. At the same time, by direct generalization of the case re = 1 we easily 
find a family of such modules for any equation of the form (|6|) . 

Proposition 10. If the coefficient of dt in any vector field from an n-dimensional involutive 
module Q is zero then Q is a singular module for the differential function L = Uf — H(t,x, ui r ^ x \). 
The co-order of singularity of Q equals one, sco£ Q = 1. 

Proof. Suppose that t s = 0. Then rank(£ s *) = n and up to changing basis of Q we can set 



All the derivatives of u with respect to x from order 1 up to order r can be expressed, on the 
manifold Q(r), via t, x and u: 



Here and in what follows 1 ^ \a\ ^ r and ao = 0. In view of the module Q being involutive, 
the vector fields Q s commute. (Moreover, we have rj s = — for some smooth function 
$ = $>(x,u) with $ u ^ 0.) Therefore, the above representation for u a is well defined since it 
does not depend on the ordering of operators in the right hand side expression. Using these 
expressions for excluding the appropriate derivatives u a from L, we get the differential function 
L = ut — H(t,x,u), where H = H(t,x,u) is the function obtained from H by the substitution 
of h a instead of u a . The order of L equals 1. Hence the vector field Q is singular for the 
differential function L, and its singularity co-order equals 1. □ 

Corollary 6. The module M = (d\, . . . , d n , d u ) is a first co-order strictly meta-singular (n+1)- 
dimensional module for any (1+n)- dimensional evolution equation. 

Proof. A p-dimensional involutive submodule Q of M satisfies the rank condition if and only 
if it is spanned by the vector fields d s — ($ S /$ U )5 U for some smooth function = <&(x,u). It 
follows from Proposition [TU] that the submodule Q is a first co-order strictly singular module for 
any evolution equation of the form Q. Varying the parameter- function we obtain a family 
of such submodules parameterized by an arbitrary function of all independent and dependent 
variables. □ 




(6) 



Q s = T s (t,x,u)d t + £ sl (t,x,u)di + r} s (t,x,u)di 



Q s = d s + f] s (t,x,u)d r , 



(7) 



Ua = h a (t, x, u) := (di + f] 



d u ) ai ■■■(dn + n n d u ) a -u. 
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The vector fields d s and d u generating the meta-singular module M commute and the differ- 
ential function L contains only first-order differentiation with respect to t (namely, in the form 
of the derivative Ut). This perfectly agrees with Theorem [TJ 

We denote by TZq(C) the set of reduction modules of £, belonging to M. Consider a module Q 
from TZq (£) and choose its basis consisting of operators Q s of the form (J7J). The system DEq(C) 
of determining equations for coefficients rj s is naturally partitioned into two subsystems. The 
first subsystem 

V S s> + V S 'r] s u = Vt' + rfvt (8) 

follows from the condition that the module Q is involutive and hence the basis elements Q s 
commute. The second subsystem 

r)l + Hr) 8 u = H s + V s H u (9) 

is a consequence of the conditional invariance criterion applied to the equation C and Q. The 
function H = H(t, x, u) is defined in the proof of Proposition [TUJ It coincides with H on 
the manifold Q (r) . The total number of equations in the joint system DEo(£) of © and © 
is n(n + l)/2 and thereby greater than the number of the unknown functions rj s if n > 1. 
Hence the system DEo(£) looks strongly overdetermined in the multidimensional case. In fact, 
the equations of the subsystems agree well with each other. The subsystem ([8j) implies the 
representation n s = — of n s via a single arbitrary function <i> = x,u). Substituting 
this representation into the subsystem ([9]), we obtain a system of p partial differential equations 
in the single function <3?, that is equivalent to a single equation in <3?. Indeed, we have 

*■ + ^-rf-<*= fl ?-£ fl ^ + (£) 1 . + (£). fl, 



i.e., this system is equivalent to the equation <I>t+$ u i?* = x(t, Here x is an arbitrary smooth 
function of its arguments and coincides with H under the substitution rj s = — $ S /$ M . The 
expression for involves derivatives of $ up to order r + 1. The function $ associated with a 
fixed module Q is defined up to the transformation = 9(t, $). As fj s = — = — = 
rj s , the functions H and H coincide. At the same time, if we choose 6 satisfying the equation 
Qt + = then $ t + $ U H® = 0. Therefore, up to the equivalence on the set of functions 
parameterizing singular modules we can assume that the function $ is a solution of the equation 

Recalling the above arguments, we derive the following result. 

Proposition 11. The system DEo(£) of the determining equations ([8]) and ([9]) is reduced by 
the composition of the nonlocal substitution n s = — $ s /$ n , where $ is a function of(t,x,u), and 
the hodograph transformation 

the new independent variables: t = t, Xi = Xj, k = 
the new dependent variable: u = u 

to the initial equation C in the function u = u(t,x, >c) with x playing the role of a parameter. 

Proof. The subsystem ([8]) implies the existence of a function <3? = <&(i, x, u) with 3> u ^ such 
that n s = — $ s /$ u . The entire system DEq(C) is reduced, up to equivalence on the set tra- 
versed by the parametric function $ from the above representation for n s , to the single equation 
&t + &uH® = 0. The last equation is mapped by the hodograph transformation to the ini- 
tial equation C for the function u = u(t,x,x). This directly follows from the definition of 
the function iJ* and the rule for calculating derivatives under the hodograph transformation, 
Ux, = -$i/$u, etc. □ 
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Proposition 12. For any equation C of the form ([6]) ; the following statements are equivalent: 

1 ) The span Q = (Q s = d s + r] s (t, x, u)d u , s = 1, . . . , n) is a reduction module of C 

2) The module Q = {Qo, Qi, . . . , Q n ), where Qo = dt + Hd u and the function H = H{t, x, u) 
coincides with H on the manifold Q(r), is involutive. 

3) There exists a vector field Qo = dt+rp(t, x, u)d u such that the module Q = {Qo, Qi, ■ ■ ■ , Q n ) 
is ultra-singular for C. 

Moreover, under these equivalent conditions the coefficient rp is uniquely determined as rp = 
H, i.e., we necessarily have Qo = Qo- 

Proof. Both the first and second statement are equivalent to the fact that the coefficients ij s 
satisfy the system DEo(£) consisting of the equations (|SJ) and ([9|). The module Q is ultra- 
singular for C if and only if rp = H, i.e. the module Q coincides with Q, and this module is 
involutive. □ 

Theorem 4. Up to the equivalence of solution families, for any equation of the form (0) 
there exists a one-to-one correspondence between one-parametric families of its solutions and 
n-dimensional reduction modules formed by vector fields with zero coefficient of dt- Namely, 
each module of the above kind corresponds to the family of solutions which are invariant with 
respect to this module. The problems of the construction of all one-parametric solution families 
of equation ([6]) and the exhaustive description of its n-dimensional reduction modules formed by 
vector fields with vanishing coefficient of dt are completely equivalent. 

Proof. Let C be an equation from class (|6|) and Q = (Qi, . . . , Q n ) € TZq(C), i.e., Q s = d s + n s d u , 
where the coefficients n s = ij s (t,x,u) satisfy the system DEo(£). An ansatz constructed with Q 
has the form u = f(t,x,ip(ui)), where / = f(t,x,ip) is a given function, f v ^ 0, ip = tp(u) 
is the new unknown function and uj = t is the invariant independent variable. This ansatz 
reduces £ to a first-order ordinary differential equation C' in ip, solvable with respect to <p u . 
The general solution of the reduced equation C can be represented in the form <p = ip(uj,z<), 
where (p H ^ and x is an arbitrary constant. The form of the general solution is defined up 
to a transformation ok = C( x ) °f the parameter x. Substituting this solution into the ansatz 
results in the one-parametric family J- of solutions u = f(t, x, x) of C with / = f(t, x, ip(t, x)), 
and any Q-invariant solution of C belongs to this family. Expressing the parameter x from the 
equality u = f(t,x, x), we obtain that x = $(t,x,u), where <3? u / 0. Then n s = u s = — <fr s /& u 
for any u £ T, i.e., for any admissible value of (t,x, x). This implies that n s = — & S /Q U for any 
admissible value of (t,x,u). 

The proof of the converse assertion is similar to that of Theorem[3j Consider a one-parametric 
family T = {u = f(t, x, x)} of solutions of C. The derivative f H is nonzero since the parameter 
x is essential. We express x from the equality u = f(t,x,x): x = <&(t,x,u) for some function 
$ = $>(t,x,u) with <I> U ^ 0. The span Q = (Qi, ■ ■ ■ ,Q n ), where Q s = d s + rj s d u and the 
coefficients n s = rj s (t,x,u) is defined by n s = — & s /<& u , is an n-dimensional involutive module. 
For any u G T we have Q[u] = 0. The ansatz u = f(t, x, <p(u})), where u = t, associated with Q, 
reduces C to the equation ip u = 0. Therefore (see [32]), Q G TZq(C). 

One-parametric families J- = {u = f(t, x, x)} and T = {u = f(t, x, x)} of solutions of C are 
equivalent if and only if the associated functions $ = $(t, x, u) and $ = x,u) satisfy the 
condition = $ S /<I> U . Therefore, equivalent one-parametric families of solutions correspond 

to the same module Q from TZq(£) and, conversely, any two one-parametric families of Q- 
invariant solutions are equivalent. □ 

Remark 3. The triviality of the above ansatz and the reduced equation results from the above 
special representation for the solutions of the determining equation. In this approach difficulties 
in the construction of ansatzes and the integration of the reduced equations are replaced by 
difficulties in obtaining the representation for the coefficients of the reduction modules. 
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Remark 4. In fact, Theorem [T] is a consequence of Theorem [3] and Proposition [TTJ This 
observation provides a uniform background for no-go results on reduction modules. It suffices to 
note that any Q-invariant solution of £ is Q-invariant, where the modules Q and Q are defined in 
Proposition [TTJ We have given a direct proof as it leads to a deeper understanding of reduction 
of evolution equations to first-order ordinary differential equations. 

8 Reduction modules of singularity co-order one 

Taking the previous example on evolution equations as a model case, we now proceed to study- 
ing (n— l)-dimensional co-order one singular reduction modules of general partial differential 
equations in one dependent and n independent variables. In the course of considering single 
modules of this kind it is only possible to derive an assertion similar to Proposition [12j 

Proposition 13. Let £: L[u] = be a partial differential equation in one dependent andn inde- 
pendent variables, let Q be an (n— 1)- dimensional involutive module of vector fields defined in the 
space of (x,u) and satisfying the rank condition, as well as wscocQ = 1- We also suppose that 
a first-order differential function L associated with L on the manifold Q (r) up to nonvanishing 
multiplier, where r = ordL, is of maximal rank with respect to the unique first-order derivative 
of u appearing in L. Then Q is a reduction module of £ if and only if there exists a (unique) 
n-dimensional module Q which is ultra-singular for £ and contains Q. 

Proof. Under these assumptions, there exists a vector field Qq such that the equation L = 
is equivalent to the equation Qo[u] = 0, where Qo[u] is the characteristic of Qq. Consider the 
module Q spanned by Q and Qq. It is n-dimensional and satisfies the rank condition. Lemma [TJ 
implies that Q is a reduction module of C if and only if it is a reduction module of the equation 
L = 0. The last condition is equivalent to the involution of the module Q, and then the 
equation C is an identity on the well-defined manifold Q(r), i-e., the module Q is ultra-singular 
for £. The conditions that the n-dimensional module Q is ultra-singular for £ and contains Q 
uniquely determine Q. □ 

Remark 5. The assumption that the differential function L is of maximal rank with respect to 
the unique first-order derivative of u appearing in L can be replaced by the weaker supposition 
that the equation L = is equivalent to the equation L = 0, where the differential function L 
satisfies the above condition of maximal rank. 

In order to have a richer theory, we shall consider families of (n— l)-dimensional co-order 
one singular reduction modules parameterized by arbitrary functions. Thus let £: L = be a 
differential equation, with L = L[u] a differential function of order r > 1. We assume that the 
function L admits a first co-order n-dimensional meta-singular module M of vector fields. (Due 
to Corollary HI we may restrict ourselves to considering only strongly singular modules of vector 
fields for differential equations.) As the special case of two independent variables when meta- 
singular modules may be non-involutive was considered in [15] . in what follows we additionally 
assume that the module M is involutive. 

Up to a change of variables we may suppose that the module M contains a first co-order 
singular family 93T = {Q*} of (n— l)-dimensional modules in reduced form parameterized by 
an arbitrary smooth function $ of (x,u), i.e., = (d s + rj s d u ), where rj s = — & S /Q U and the 
index s runs from 1 to p = n — 1. We use the representation for r] s from the very beginning. 

By Theorem [TJ the differential function L can be written in the form L = L{x, ^ r ,i,n-i)j 
where J7 r ,i,n— l = { u a, ot n ^ 1, \a\ ^ r), L Ua ^ for some u a with a n = 1. In this case the 
restriction of L to Q%, coincides with the restriction, to the same manifold Q%\, of the function 
L = L(x,Q, rt i tn -i), where 

O r , 1 , n _ 1 = (££"(Qfr ■■■{Qt-i) an - 1 u,a n < 1, |a| < r). 



23 



Consequently, the form of L is determined by the form of L and a chosen value of the parameter- 
function Depending on the value of <£, the differential function may either identically 
vanish or be of order or 1. Thus either the module Q* is ultra-singular or sco^ = or 
scolQ* = 1- We analyze each of these cases separately. In addition, we suppose that the 
function L* is of maximal rank with respect to u (resp. u n ) if sco^ = (resp. sco^ = 1). 

The condition L* = 0, where u and u n are considered as independent variables, determines 
those values of <3> for which the module is ultra-singular for C. We split this condition 
with respect to u n , thereby obtaining a system <S_i of partial differential equations in <I> of 
orders not greater than r, which may be incompatible in the general case. Incompatibility here 
amounts to the family 9JT containing no ultra-singular modules. E.g., evolution equations of 
orders greater than 1 have no ultra-singular modules generated by vector fields of the general 
form ([7]), see Section [71 That <I> satisfies the ultra-singularity condition S-\ guarantees that 
€ 7£ n-1 (>C) and the family of Q*-invariant solutions of C is parameterized by an arbitrary 
function of the single Q^- invariant variable x n . 

Under the assumption sco^Q* = it follows that the parameter- function $ satisfies the 
condition L* = with u and u n viewed as independent variables, which is weaker than the ultra- 
singularity condition. In this case the corresponding system Sq of partial differential equations 
in $ of orders not greater than r, obtained by splitting the zero co-order singularity condition 
with respect to u n , is more likely to be compatible than <S_i. 

Next we mention certain sufficient conditions for the compatibility of Sq. If L Un = then 
the system So is compatible since it is satisfied by any with ($ S /$ U ) M = 0. Up to equivalence 
of functions parameterizing modules from the family we can assume that $ = u — C{x) and 
r] s = ( s , i.e., Q u ~^ = (d s + ( s d u ). In other words, we have sco^ Q" - ^ ^ for any £ = ((x). 

Assuming L u = 0, the condition L u ~^ = under the assumption £ = £( x ) implies only a 
single partial differential equation with respect to £. Any of its solutions is a solution of <S_i 
and hence the corresponding module Q u ~^ is ultra-singular for L. 

Otherwise sco^ Q u ~^ = and we can resolve the equation C with respect to u as a variable of 
the underlying jet space. In other words, we represent this equation in the form u = K[u]. Here 
K[u] is a differential function which depends at most on the variables x and u a with < |or| ^ r, 
a n ^ 1 and, if a n = 1, \a\ > 1. The equation L u ~^ = is obtained via replacing u s by (, s (x) 
in C Therefore, it can be represented in the form u = G^{x), where G^{x) = K[v\\ u= ^i x y Thus, 
the expression for the function involves derivatives of the parameter- function £ = up 
to order r. The conditional invariance of the equation C with respect to the module Q u ~^ is 
equivalent to the compatibility of the system u = G^, u s = ( s with respect to u and hence 
leads to the system of n — 1 partial differential equations ( s = G| with respect to £. Here 
G^ is interpreted as a differential function of £ = C( x )- Then d is the total derivative of G^ 
with respect to x s . As the function-parameter £ is defined for a fixed module up to a constant 
summand, the system is equivalent to the single partial differential equation £ = G">. This means 
that Q u ~^ is a reduction module of the equation C if and only if up to shifts of the dependent 
variable the parameter-function £ is a solution of the same equation. The ansatz constructed 
with Q u ~^ can be taken in the form u = tp(uj) + C( x )} where ip = ip{oj) is the new unknown 
function and oj = x n is the invariant independent variable. It reduces the initial equation C to 
the trivial algebraic equation (p = 0, i.e., the function u = Ci x ) 1S the unique (^"''-invariant 
solution of C Conversely, we fix a solution u = £(x) of the equation C Then £ = G ( '(x) and 
hence £ s = Gs, which is equivalent to the conditional invariance criterion for the case of the 
module Q" - '' = (d s + Csd u ) and the equation C, i.e., Q u ^^ is a reduction module of C, and 
Q u = 0. The solution u = C( x ) is invariant with respect to Q u ~^ by construction. Thus we 
obtain: 

Theorem 5. Suppose that an equation C: L = possesses a first co-order singular family SOT = 
{Q*} of (n—1) -dimensional modules in reduced form = (d s — (3> s /§ u )d u ) parameterized by an 
arbitrary smooth function $ of(x,u), i.e., the left hand side L of this equation is represented in 



24 



the form L = L(x,Q rj i :n -i), where JV,i,n-i = { u a, ct n ^ 1, [a| ^ r), L Ua ^ /or some u a with 
a n = 1, and additionally L Un = and L u ^ 0. T/ien there exists a one-to-one correspondence 
between solutions of C and reduction modules from 5DT with <3? = u — C( x )- Namely, any such 
module is of singularity co-order and corresponds to the unique solution which is invariant with 
respect to this module. The problems of solving the equation L and the exhaustive description of 
its reduction modules of the above form are completely equivalent. 

Next we turn to analyzing the regular values of $ for which the singularity co-order of Q* 
coincides with the singularity co-order of the whole family 9Jt (and equals 1). If sco^Q* = 1, 
the parameter- function $ satisfies the regularity condition L® n ^ 0. Thus the equation = 
which is equivalent to C on the manifold Qf r) can be solved with respect to u n : u n = G®(x,u). 
Here the expression for the function G* depends on derivatives of the parameter-function $ 
up to order r. The conditional invariance criterion, when applied to the equation C and the 
operator G;*, gives the system 

< + <G* = G? + r, s G* (10) 

with respect to the function $ (recall that r/ s = — <3? S /<I> M ). This system coincides with the 
formal compatibility condition of the equations u n = G* and u s = r] s with respect to u. As 
the function G* can be also expressed directly in terms of the coefficients rj s , the system (|10p 
supplemented with the involution condition if s , + r] s 77* = r/f + r/ s r/* can be interpreted as a 
system with respect to r] s , cf. Section [71 

In fact, the system (|10p is equivalent to a single equation in <£. Indeed, 

<* + i-Ct - < - < * - Gf - + (1)^ + (£) <G . 

= ¥:( 8 --l: 8 ") ($ " +$ " o4)=0 - 

i.e., system (fit)]) is equivalent to the equation <!>„ + $„G $ = x(x n ,$), where x is an arbitrary 
smooth function of its arguments. The function <3? associated with a fixed module G;* is defined 
up to the transformation $ = 9(x n , 4>). As 77 s = — $ s /$ u = — $ s /$ u = 77 s , the functions G* 
and G* coincide. At the same time, if we choose 9 satisfying the equation 6 n + = then 
3*n + ^*nG* = 0. Therefore, up to the equivalence on the set of functions parameterizing singular 
modules we can assume that the function $ is a solution of the equation <3> n + $ M G* = 0. 

The order of each equation from system ()10p with respect to $ equals r+1 and hence is greater 
than the order of the system <So- Under certain smoothness assumptions (e.g., analyticity) this 
implies that the equation (jlpp has solutions which are not solutions of Sq. Consequently, the 
equation C necessarily possesses first co-order singular reduction modules which belong to 9Jt. 

Proposition 14. Suppose that any (n—1) -dimensional involutive module of the reduced form 
Q = (d s + rj s d u ) is a first co- order singular module for an equation C Then the determining 
system for values of rj s associated with reduction modules of C is reduced by the composition of 
the nonlocal substitution r] s = — § S /Q U , where § is a smooth function of (x,u) with $ u 7^ 0, and 
the hodograph transformation 

the new independent variables: Xi = x%, k = <3?, 
the new dependent variable: u = u 

to the initial equation C for the function u = u(x, x) with k playing the role of a parameter. 

Proof. The possibility of representing rj s in the form §> s /$ u with some function $ = <&(x,u) 
follows from the involution condition r/f, +rj s rj^ = +rj s r]f l for the module Q. In this represen- 
tation, the system of determining equations for values of rj s associated with reduction modules 
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of C has the form (fit)]) and is equivalent, up to equivalence on the set traversed by the parametric 
function to the single equation <J> n + Q U G® = 0. In view of the definition of the function G* 
and the expressions for derivatives under the hodograph transformation, % j = — $j/$ u , etc., 
the hodograph transformation maps the equation $ n + $ U G* = into to the initial equation C 



The above connection between the initial equation C and the determining equation (jlOp can be 
stated as a relation between one-parametric families of solutions and (n— l)-dimensional first co- 
order singular reduction modules. The results of Section [5] imply that for each (n— l)-dimensional 
first co-order singular reduction module Q of the equation C there exists a one-parametric family 
of Q-invariant solutions of C. If the equation C admits an n-dimensional meta-singular module 
of singularity co-order one, the converse statement is true as well. It is convenient to prove this 
statement without transforming the meta-singular module to the reduced form. 

Theorem 6. Suppose that an equation C possesses an n-dimensional first co-order meta-singular 
module M . Then for any one-parametric family J- of solutions of C there exists an (n— 1)- 
dimensional involutive submodule Q of M that is a reduction module of C and each solution 
from T is Q-invariant. 

Proof. Let J- = {u = f(x, x)} be a one-parametric family of solutions of C. Here, f K is nonzero 
since the parameter x is essential. Prom u = f(x,x) we conclude that n = &(x,u) with some 
function $ = <b(x,u) with $ n / 0. 

Let {Qo, ■ ■ ■ ,Q n } be a basis of M. Suppose that there exists a 6 {0, ...,n} such that 
Q a § ^ 0. Up to permutation of basis elements we can assume that Qo& 7^ 0. Then we set 
Qs = Qs~ (Qs$)/(Qo$)Qo- If Qa$ = for all a G {0, . . . , n}, we set Q s = Q s . 

Consider the submodule Q generated by the vector fields Q s . This submodule is involutive 
and of dimension n — 1. Hence scolQ* ^1. It is also obvious that Q S Q = 0. As /j = 
— (&i/& u )\ u= f, this means that any solution from the family T is Q-invariant. The case sco^ Q = 
is impossible as otherwise the equation C could not have a one-parametric family of Q-invariant 
solutions. Therefore either Q is an ultra-singular module for L or sco/, Q = 1. Any ultra-singular 
module for L is a reduction module of C If sco^ Q = 1 then Q is a reduction module of C by 
Proposition [HJ □ 

For the correspondence between one-parametric families of solutions and (n— l)-dimensional 
reduction modules to be one-to-one, the related meta-singular module should satisfy additional 
restrictions. 

Theorem 7. Suppose that an equation C possesses an n-dimensional first co-order meta-singular 
module M and all (n— 1)- dimensional submodules of M as well as the entire module M are not 
ultra- singular for C Then up to the equivalence of solution families there exists a bijection 
between one-parametric families of solutions of £ and its (n— 1)- dimensional reduction modules 
contained in M . Namely, each module of this kind corresponds to the family of solutions which 
are invariant with respect to it. 

In other words, the problems of the construction of all one-parametric solution families of 
the equation C and the exhaustive description of its reduction modules of the above form are 
completely equivalent. 

Proof. If Q is an (n— l)-dimensional reduction module of C contained in M, we have sco£ Q = 1. 
Therefore, in view of Proposition [8] the equation £ possesses a one-parametric family T of Q- 
invariant solutions, and any Q-invariant solution of C belongs to this family. Each one-parametric 
family of Q-invariant solutions of £. is obtained from T by re-parameterizing. 

Conversely, consider a one-parametric family T = {u = f(x, x)} of solutions of L. Theorem^] 
implies that there exists Q € TZ n ~ 1 (C) such that Q is a submodule of M and any solution 
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from J- is Q-invariant. Let us prove that the module Q is unique. Suppose that there exists an 
(n— l)-dimensional involutive submodule Q of M different from Q and any solution from T is 
Q-invariant. This implies that the family T consists of solutions invariant with respect to the 
entire module M. 

Therefore, the module M satisfies the rank condition. To show this, we fix a basis 
{Qo, ■ ■ ■ ,Q P } of M, where Q a = £ CTl (x,u)<9j + rj a (x, u)d u , a = 0, ...,ra. The function <I> de- 
fined in the proof of Theorem [6] is an invariant of all Q a , i.e., Q a Q = 0. As <1> U ^ 0, we have 
that rank(^ (TJ ) = n. 

As M is an n-dimensional involutive module satisfying the rank condition and the family J- 
formed by M-invariant solutions of the equation £ is one-parametric, in view of Proposition [7J 
the module M is ultra-singular for C, thereby contradicting an assumption of the theorem. This 
means that the module Q is unique. □ 

The bijection discussed in Theorem [7] is broken in the presence of (n— l)-dimensional ultra- 
singular submodules or ultra-singularity of the entire meta-singular module. Indeed, if an (n— 1)- 
dimensional involutive module Q is ultra-singular for the equation C, the family of Q-invariant 
solutions of C is parameterized by an arbitrary function of a single argument, cf. Proposition [7J 
If the entire n-dimensional module M is ultra-singular for the equation C, this equation possesses 
a one-parametric family T of M-invariant solutions. Then any solution from the family T is 
invariant with respect to each (n— l)-dimensional involutive submodule of M. 



9 Singular modules for quasi-linear second-order PDEs 

For some classes of differential equations it is possible to exhaustively describe the associated 
singular modules. We study certain quasi-linear second-order PDEs from this point of view. 
It is natural to distinguish elliptic, evolution and generalized wave equations. In all cases, Q 
is an involutive module of vector fields defined in the corresponding space of independent and 
dependent variables and satisfies the rank condition, and the dimension of Q is less than the 
number of independent variables. 

Elliptic equations. Consider an equation C for the single unknown function u of the indepen- 
dent variables x = [x\, . . . , x n ), having the general form 

L[u] := a lj (x,u {Vl )uij + b(x,u (1) ) = 0, (11) 

where the coefficients a 13 and b are defined on the same domain O of the first-order jet space, 
a iJ = a 31 and the matrix- function {a 13 ) is positive definite in each point of Vt. We will prove that 
the equation L possesses no singular modules of dimensions less than n. 

Denote dim Q by p, < p < n. Up to permutation of x we can locally choose a basis of Q 
which consists of vector fields of the form Q s = d s + £ st (x, u)d L + f) s {x, u)d u , cf. Section^ Here 
and in what follows the index i runs from p + 1 to n. Then any derivative of u of order one or 
two is expressed, on the manifold Q ( 2), via derivatives of u with respect to x p+ \, . . . , x n only 
and the coefficients £ st and 77 s , see equation ([1]). As only second-order terms in the expressions 
for second-order derivatives of u are essential here, we use the representations 

u u = -i SL 'u u > + R SL (x,u (1) ), u ss > = C£ sU 'u u , + R ss '(x,u (1) ), 

where R sl denote the terms without second-order derivatives. Substituting the above expressions 
for the derivatives u s i into L[u], we obtain the differential function L[u] := a u u u > + b, which is 
associated with L[u] on the manifold Q(i)- Here a" = a u — a st ^ st +a ss L and b = b+a SL R SL + 
a ss R ss . In other words, for each fixed 1 and 1! = 1 the coefficient a u coincides with the value of 
the quadratic form whose matrix is (a 13 ) at the tuple (z , . . . ,z n ), where z s = — z l = 1 and the 
other z's equal zero. As the matrix {a 13 ) is positive definite, the coefficient (a 13 ) is nonvanishing. 
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This implies that the differential function L[u] cannot coincide on the manifold Q,m with a 
differential function of order less than two, not even up to a nonvanishing multiplier. Therefore, 
wsco£ Q = sco£ Q = 2. 

Evolution equations. Similarly to Section [JJ for this and the next class of equations we 
set the number of independent variables to n + 1 instead of n and additionally single out the 
variable t = xq, i.e., the unknown function u depends on the variables t and x = (x\, . . . ,x n ). 
The general form of a quasi-linear second-order evolution equation C we intend to study is 

Ut = H[u] := a lj (t,x,u {1:X) )uij +b(t,x,u (ltX) ), (12) 

where uq_ >X ) = (u,ui, . . . ,u n ), the coefficients a ij and b are defined on the same domain Q as 
traversed by their arguments, and the matrix-function {a 13 ) is symmetric and positive definite in 
each point of fi. Up to permutation of the variables X(, any involutive module Q of vector fields 
defined in the space of (t,x,u) and satisfying the rank condition, where dimQ < n + 1, can be 
locally assumed to be spanned either by the vector fields &t + i 0L d L + rpd u , d s + £ SL d L + i] s d u with 
p = dim Q — 1 or by the vector fields d s + ^ SL d L +rj s d u with p = dimQ. Here £ 0t , rf , £ st and r] s are 
smooth functions of (t,x,u). We can show similarly as for elliptic equations that in the second 
case with p = n and only in this case the module Q is singular for C, and wsco£ Q = sco£ Q = 1. 
This is the case that has been studied in Section [7J 

Generalized wave equations. Consider the equation £: uu = H[u] for the single unknown 
function u of the independent variables t and x = (x\, . . . ,x n ), where the differential function 
H = H[u] is defined analogously to (fT2|) but the coefficients a 13 and b may additionally depend 
on lit- We partition the set of appropriate involutive modules in a way different from that used 
for evolution equations. Up to permutation of the variables x, any involutive module Q of vector 
fields defined in the space of (t,x,u) and satisfying the rank condition, where dimQ < n + 1, 
can be locally assumed to be spanned either by the vector fields dt + t]°d u , d s + £, SL d L + rj s d u with 
p = dimQ — 1 or by the vector fields d s + T s dt + £ st <9 t + r] s d u with p = dimQ. Here 77°, r s , £ st 
and 77 s are smooth functions of (t,x,u). Only in the second case with p = n the module Q may 
be singular for C. We consider this case in detail. 

As p = n, the basis elements of Q take the form Q s = d s +T s dt+r] s d u . On the manifold Q ( 2) we 
have that u s = r/ s — £ s ut, u ss i = £ s £ s u# + R ss , where R ss denotes the corresponding collection 
of the terms depending at most on t, x, u, ut and uu- Substituting the above expressions for the 
derivatives of u with respect to x into the differential function L[u] = uu — H[u], we obtain the 
differential function L[u] := (1 — a 1 - 3 t 1 t 3 )uu + b(t, x,u,ut), which is associated with L[u] on the 
manifold Q ( 2) and depends at most on t, x, u and ut- Here the coefficient a 13 = a 13 (t,x,u,ut) 
is obtained from a 13 by the substitution u s = r] s — £ s ut and the precise form of the coefficient 
b = b(t,x,u,ut) is not essential. The module Q is singular for £ if and only if ordL[n] < 2, i.e., 
the coefficients t s = r s (t,x,u) satisfy the equation o 13 t 1 t 3 = 1. If some of the coefficients a 13 
depend on Ut, this equation should be split with respect to this derivative and hence it may be 
inconsistent. 

As the module Q is involutive and, therefore, the vector fields Q s commute, the coefficients 
t s and r] s jointly satisfy the system Q s t s = Q s t s , Q Sr q s = Q s rj s , i.e., 



In view of the Frobenius theorem this implies that the system = with respect to the 
unknown function <3? = f>(t, x, u) admits solutions I = 1,2, such that 3>£$^ — $4$? 7^ 0. 
Solving the pair of the equations & s + T s & t + r] s & u = for each fixed s as a system of linear 
algebraic equations with respect to (T s ,r] s ), we derive the representation 



rf + T s rf + rj s r^ = t s s , + T s 'r t s + 

+ r s ri s t ' + rf = 7$ + T s ' V s t + rf 'r/*. 



(I>1(T)2 _ &2fo 
_ s u _ s _ 



(13) 
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Conversely, it can be checked by direct calculation that for arbitrary functions $ 1 and <E> 2 with 
&}<&u ~ &u^t 7^ the module Q spanned by the vector fields Q s = d s + r s d t + r] s d u with 
the coefficients defined by (|13p is involutive. Up to functional dependence of pairs of functions 

and <£ 2 with <I> 2 — 7^ 0, there exists a bijection between such pairs and involutive 

modules spanned by n vector fields of the form Q s = d s + r s dt + rj s d u . 

If the coefficients a lJ do not depend on ut (this is the case if, e.g., o lJ = a l i (t,x,u)), the 
substitution of the expressions (|13j) for t s into the equation a^rV- 7 = 1 gives a single equation 
in the two unknown functions and <J> 2 . Each solution of this equation is associated with 
a singular module Q for the equation C. In general, in the multidimensional case n > 1, 
the coefficients t s and rj s of the corresponding basis elements Q s are coupled in a nonlocal 
and nonlinear way. This is why it is impossible to exhaustively describe singular modules of 
multidimensional nonlinear wave equations within the framework of meta-singular modules. At 
the same time, under additional conditions for the coefficients a tJ nonlinear wave equations 
possess families of singular modules which fit well into the above framework. 

Thus, let the coefficients a lJ only depend on t and x. We look for singular modules of jC for 
which the corresponding coefficients r s also do not depend on u. It is then sufficient to consider 
pairs of functions and <3? 2 with <3?^ = and, therefore, <E> 2 7^ 0. Then the expressions (|13j) 
for t s are reduced to t s = —$>1/<$>1. The equation a !J rV J = 1 is equivalent to the equation 
((J>!) 2 = a l i$>j&j, which is the eiconal equation associated with the wave equation C We fix a 
solution ^ = ^(t,x) of the eiconal equation with vl/ t 7^ and consider the module M* spanned 
by the vector fields ^td s — ^ s dt, s = 1, . . . , n, and d u . The module M* is meta-singular for the 
equation C, and sco£ ^ 1 for each ^. Hence, the results of Section[8]are relevant in this case. 
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